Theory of Pendular Rings Revisited 



Boris Y. Rubinstein 1 and Leonid G. Fel 2 * 

^towers Institute for Medical Research, 1000 E 50th St, Kansas City, MO 641 10, USA 
2 Department of Civil and Environmental Engineering, 
Technion - Israel Institute of Technology, Haifa, 32000, Israel 

August 1, 2012 



Abstract 

We present the theory of liquid bridges between two axisymmetric solids, sphere and plane, with 
prescribed contact angles in a general setup, when the solids are non-touching, touching or intersecting, 
We give a detailed derivation of expressions for curvature, volume and surface area of pendular ring as 
functions of the filling angle ip for all available types of menisci: catenoid Cat, sphere Sph, cylinder 
Cyl, nodoid Nod and unduloid U nd (the meridional profile of the latter may have inflection points). 

The Young-Laplace equation with boundary conditions can be viewed as a nonlinear eigenvalue 
problem. Its unduloid solutions, menisci shapes z^(r) and their curvatures H^(ip), exhibit a discrete 
spectrum and are enumerated by two indices: the number n of inflection points on the meniscus merid- 
ional profile M. and the convexity index s = ± 1 determined by the shape of a segment of A4 contacting 
the solid sphere: the shape is either convex, s = 1, or concave, s = — 1. 

For the fixed contact angles the set of the functions H^ (ijj) behaves in such a way that in the plane 
{ip, H} there exists a bounded domain where H^(ip) do not exist for any distance between solids. The 
curves H^(ip) may be tangent to the boundary of domain which is a smooth closed curve. This topolog- 
ical representation allows to classify possible curves and introduce a saddle point notion. We observe 
several types of saddle points, and give their classification. 
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1 Introduction 



The problem of pendular ring (PR) arises when a small amount of fluid forms an axisymmetric liquid bridge 
with interface (meniscus) between two axisymmetric solids. This problem includes a computation of liquid 
volume V, surface area S and surface curvature H and was one of gems in mathematical physics of the 19th 
century. In the last decade the PR problem became again an area of active research due to investigations 
on stability of the PR shapes, and its importance has grown for various applications in soil engineering, 
physics of porous media, etc. 

The history of the problem dates back to 1841 when Delaunay [1] classified all non-trivial surfaces of 
revolution with constant mean curvature in M 3 by solving the Young-Laplace (YL) equation and showed 
that they are obtained by tracing a focus of a conic section when rolled on a line, and revolving the resulting 
curve around the axis of symmetry. These are cylinder (Cyl), sphere (Sph), catenoid (Cat), nodoid (Nod) 
and unduloid (Und). The two last of them are defined through the elliptic integrals and may appear of two 
kinds, concave (-) and convex (+), depending on constant sign of the meridional profile M. curvature. One 
more type of meniscus, an inflectional unduloid, appears when meridional section Ai curvature changes its 
sign along the meniscus. 

In 1864 Plateau [6] applied this classification to analyze the figures of equilibrium of a liquid mass, and 
was the first who discovered Q a standard sequence of meniscus evolution observed with increase of the 
liquid volume in absence of gravity. According to [5 ] the Plateau sequence reads: 

Nod~ -> Cat -> Und -> Und^ -»■ Und+ -»■ Sph -> Nod+ , (1.1) 

where subscripts denote the number of inflection points on the meniscus meridional section M. . Even so, 
the actual algorithm for solution of the PR problem leads to the eigenvalue problem for mean curvature H 
that requires extensive and accurate computation of the elliptic integrals and was not available before the 
computer era has been started. A complete review on different methods used to find actual solutions of the 
YL equation or the equivalent variational problem (Howe in 1887 and Fisher [2] in 1926) throughout 
the last century can be found in Q. 

In 1966 Melrose (H gave a detailed analysis of the Nod - meniscus and derived the formulas for V, S 
and H in the case of two touching spheres of equal radii. Orr, Scriven and Rivas in 1975 extended this result 
in seminal article [5] for the menisci of various profiles in the case of solid sphere of radius R above the 
solid plane for d = D/R > with prescribed contact angles Q\ and 6*2 on sphere and plane, respectively; 
here D denotes a distance between the sphere and the plane. In the case of touching solids, d = 0, they (21 
have performed numerical computations and verified the Plateau sequence (11.11 ) when the liquid volume is 
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increasing. 

During the past decades the work Q became classical, albeit throughout a vast number of references 
no attempt was made to extend this explicit analysis using modern computer algebra technique. What is 
more important, formulas in [5] for the non-touching solids were left without detailed analysis. The reason 
why we have noticed this fact is based on substantial difference in the behavior of the function H(ip,d) in 
two different setups: d = and d > 0. This can be seen easily in Figure Q] where we consider for small 
filling angles ip two types of menisci between touching (a) and non-touching (b) sphere and plane having 
the same set of contact angles 9\ < it/2 and 02 = vr/2. 






. 







(a) (b) 

Figure 1: (a) The shape of meniscus for d = 0, 0\ = vr/6, 62 = w/2, ip = vr/30, H a R = —69.57. (b) The 
shape of meniscus for d = 0.076, 9\ = 7r/6, 62 = vr/2, ip = vr/30, H^R = 2.14. 

In the case (a) the sphere-plane geometry approaches its wedge limit while in the case (b) it approaches 
the slab geometry. Estimate two principal radii, meridional R v and horizontal R^ in both cases for ip -C 1. 
In the first case (a) they are of different signs, i.e., R v < 0, R h > 0. Keeping in mind H = l/2(R- 1 +R' h 1 ), 
a simple trigonometry gives 

R 2 cos V\ R V 
The dependence H ~ — ip~ 2 , tp <C 1 describes the asymptotics H(tp) of the Nod - meniscus for two 
touching solids found in @. In section IC2l we justify formula (11.21) by rigorous derivation of nodoidal 
asymptotics. In the second case (b) we have another estimate, 

^-j- , RH^^j. (1.3) 

D cos 61 R 2ip 

We have found meniscus with H > that according to the Plateau classification has the concave unduloid 
type Undo which is absent in the range ip -C 1 in the sequence (ll.lb . Here another asymptotics holds, 
H ~ ip~ l . This leads to dramatic changes for the whole sequence (11.11 ) giving rise to the Cat menisci for 
two different filling angles ip, or to a single degenerated Cat meniscus, or to disappearance of both of them. 

There exists one more case which cannot be reduced to the previous ones. This is a meniscus between 
intersecting sphere and plane (d < 0) having the same set of contact angles Q\ and 62 = vr/2 and wedge 



geometry for small filling angles ip. The calculation gives two principal radii, R v and Rh, and its curvature 
as follows (we refer to the Figure \l2\ b) in section |7]), 

~ET tfl t* i \ ' ~ET - sin ^* ' RH ~ - C °nn 1+ ^\ ' ^* = arcco^ 1 + d ) • C 1 - 4 ) 

cos (0i +ip*> R 2(tp — ip*) 

The dependence H ~ — ("0 — t/j*)~ ,i[> — ip# <C 1 describes the asymptotics i? (VO of the Nod - meniscus 
for two intersecting solids and is intermediate between d 1 -2b and d 1 -3b - Thus, even a change of a single 
governing parameter d only, when other two 0i , 02 are fixed, changes drastically the evolution of menisci. 

Our analysis of solutions of the YL equation shows that the changes become more essential (non- 
uniqueness of solutions, excluded domain B in the {if}, H} plane where H(ip) does not exist, etc.) when 
we deal with the whole 3-parametric space P 3 = {0i,02,d}. From this point of view the article [5] has 
dealt with 2-parametric subspace {0i, 02, 0}. This creates an additional challenge to describe the menisci 
in different areas of P 3 , i.e., to give a complete theory. 

In this paper we present the theory of pendular rings located between two axisymmetric solids, sphere 
and plane, in a general setup, when the solids are non-touching, and touching or intersecting. We give a 
detailed derivation of expressions for curvature H(tp), volume V{ip) and surface area S(ip) as the functions 
of the filling angle for all available types of menisci including those omitted in Q. We give also an 
asymptotic analysis of these functions in the vicinity of singular points where they diverge. 

The paper is organized in eight sections and four appendices. In section [2] we give a setup of the 
problem and derive the YL equation and its solution through the elliptic integrals. The integrals introduced 
in this section are evaluated in section 12.11 and the explicit expressions for the meniscus curvature, shape, 
volume and surface area are found. We introduce also a new function a(ip) which is intimately related to 
the curvature H(tp) and becomes a main tool in analysis of an evolution of pendular rings. In section [3] we 
discuss the general curvature behavior for different types of menisci. Existence of catenoids in the menisci 
sequence for the cases of non-touching solids is considered in section 13.11 

In contrast to the case of touching solids discussed in [5 ], when for the small ifi the Nod - menisci exist, 
in the case of non-touching solids the Und^ menisci come first. This allows existence of two catenoids in 
the menisci sequences, while in some cases catenoids do not appear at all. We find the critical value of 
the distance between solids at which the catenoids merge and estimate the curvatures for different types 
of menisci. We show that the values of a can be bounded for some types of the menisci. When these 
bounds are crossed, the corresponding menisci transform one into another. Analysis of these transitions, 
their sequence and smoothness, is given in section 0] 

In section |5] we elaborate a topological approach to study different curves a n (ip) enumerated by the 
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number n of inflection points at menisci. A behavior of these curves in the plane {ip, a} is confined within 
domain A' = {0 < t/j < it — 9\, < a < 1} with embedded subdomain B C A' which is prohibited for 
a n (ip) to pass through; this makes A' not simply connected. The curves a n (tp) may be tangent to the 
subdomain's boundary which is a smooth closed curve described by symmetric transcendental function. 
This global representation allows to classify possible curves and introduce a saddle point notion in the PR 
problem. We observe several types of saddle points, their classification is presented in section [6] 

In section |7] we give a brief analysis of menisci evolution in the cases of touching and intersecting solid 
bodies which is essentially different from a general setup of non-touching bodies. Concluding remarks and 
open problems are listed in section [8] 

Four appendices are inseparable parts of the paper. Appendices lAl and iBl contain a list of formulas with 
technical details for H, V, S and the shape of meniscus of each type separately. They give an exhaustive 
description of menisci and build a basis for further investigation of basic properties of PRs like stability, 
rupture, hysteresis etc. In appendix[C]we show that the nodoid meniscus Nod + always has a local minimum 
of curvature and local maxima of the surface area and volume. We also find the asymptotic behavior of the 
nodoidal curvature in vicinity of singular point when — 2 < d < 0. Appendix 151 is completely devoted to 
elliptic integrals and their applications to computation various expressions throughout the paper. 



2 Young-Laplace Equation and its Solutions 

The problem of PR can be posed as a search of the surface of revolution characterized by a constant mean 
curvature H satisfying the YL equation valid in case of negligibly small gravity effect 

z" z' 
2H = ^ + m , (2.1) 

(1 + ^/2)3/2 r ( 1 + 2 /2)V2 

where z[r) and r are cylindrical coordinates of the meniscus. Introducing new variables x = r/R and 
y = z/R and a parameter u = s'mt (where t is an angle of the normal to meniscus with the vertical axis), 
we transform this equation into the problem for nondimensional curvature H = RH 

2H = du/dx + u/x. (2.2) 

The contact angles with the solid bodies are 6\ (with the sphere) and 02 (with the plane). The boundary 
conditions read 



ti = 0\ + ip, yi = 1 + d — cos if), x\ = sin rp, 
t 2 = ir - 9 2 , y 2 = 0. 
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(2.3) 



Here ip is the filling angle, and d = D/R is the scaled distance between the sphere and the plane. It is easy 
to show that 

dy/dx = tant . (2.4) 



The solution in parametric form reads 



1 

2H 

-f 



sin t + sV sin t + c 



sini + 



s sin t 
\J sin 2 t + c 



where we used the relation dx/dt = xs cos t/ y sin t + c and the parameter c depends on curvature 

c = 4//sin^ (H simp — sinti) . 

Here and below s = ±1; its computation will be described in section|4j formula d4.ll ). 
Introducing a parameter a, such that 

H s'mifj = as'mti, 

we rewrite relation (12.71) as 



(2.5) 
(2.6) 



(2.7) 



c = 4a(a — 1) sin 2 t±. 
Making use of the boundary conditions we find for the curvature 



(2.8) 



(2.9) 



2H^> = L 



sint + 



s sin t 



dt , ^ = d + 1 - cos ip 



(2.10) 



' j \J sin 2 t + c 

The meniscus surface area S is computed as S = 2ir J x\J\ + (dx/dy) 2 dy and is given by the integral 

2 



7T 



2tf 2 



^ , = S 



i sin t ( sin t + s y sin 2 1 + c 



*2 | sint | v sin 2 1 + c 



(it . 



(2.11) 



The volume V r = it f x 2 dy of the solid of rotation reads 



7T 



t\ sint (sint + s\/ sin 2 1 + 



dt . 



(2.12) 



'h V sin 2 t + c 

Then the volume V of liquid inside the PR can be computed by subtracting from the above expression the 
volume V ss = ir(2 — 3 cos ip + cos 3 ip) /3 of the spherical segment corresponding to the filling angle ip. 
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2.1 Integral Evaluation 

Consider evaluation of three integrals I s , J s and K s and give expressions for H, V, S in terms of c, filling 
angle ift, contact angles B\, 6> 2 and distance d from the sphere to the plane. The integral I s can be written as 
I s = I\ + sl 2 , where 

rh 

h(ti,t2) = / sint dt = — costi + cost 2 , (2.13) 
Jt 2 

r / f l1 sin 2 tdt „ „ 

/ 2 (*i,*2) = / , = (2-14) 
J*2 v sin t + c 

Making use of the elliptic integrals of the first Fit, k) and the second E(t, k) kind, respectively, we find 
(see Appendix ID. II ) 

h(ti,t 2 ) = y/c \E(ti, k) - F(ti, k) - E(t 2 , k) + F(t 2 , k)] , k 2 = — 1/c, (2.15) 

where by A(z) we denote a complex conjugation of complex function A(z). 

Consider the integral K s required for computation of the surface area in (12. lib - The value of shxt\ at 

the upper limit t\ can take both positive (for t\ < it) and negative (for t\ > n) values. In the first case the 

integral reads K s = 5(2/2 +I3) +2Ii, where 

fti c dt _ — 

h(h,t 2 )= =V~c \F(t lt k) - F(t 2 ,k)] . (2.16) 

Jt 2 v sin t + c 

In the last case the integral is broken into two parts as follows 



K s {tx,t 2 ) = K s (K,t 2 ) + K s {K,t 1 ) . (2.17) 

The integrals I2 and I3 follow another relation 

h{ti,t 2 ) = I 2 (7T,t 2 ) - / 2 (vr,ti), h(t h t 2 ) = h(ir,t 2 ) - iaMi) . (2.18) 

Using (I27T51 — l27T8T> we find for t\ > vr 

K s {h,t 2 ) = s(2/ 2 (7r,t 2 ) + 2/ 2 (vr,t 1 ) + I 3 (TT,t 2 ) + h(TT,h)) + 2/i(7r,t 2 ) + 2/ 1 (vr, t x ) 
= s(2I 2 (t u t 2 ) +h{ti,t 2 )) + 27i(ti,t 2 ) + 4(costi + 1) 
+ 2s^/c [AE(k) - 2K{k) - 2E(t!, k) + F(t u k)] . 

Finally we have 

s(2/ 2 + I 3 ) + 2h, tl<TT, 
K s (h,t 2 )=l s (2/ 2 + I 3 ) + 2/i+4(costi + l)+ (2.19) 
2 Sx /c [4E(k) - 2K(k) - 2E{t u k) + F(t h k)] , h > vr. 



Both surface area S(ip) and its derivative S'(ip) are continuous at the matching value ip = tt — 9\. 

Rewrite the integral J s required for the volume computation in (12.121) J s = 4J3 + cl\ + s{J\ + 3J2), 
where 

f* 1 sin 4 * (it ( tx . 2 r~2 , T f h ■ 3 , 

Ji = =, J 2 = sm *vsm t + cdt, J3 = sm t at. 

Jt 2 v sin 2 i + c •/ 1 2 J*2 

The last integral reads 



We find J\ = J 2 — CI2, where 
J 2 (ti,t 2 ) = 



cos 3ii — 9 cos t\ cos 3i 2 — 9 cos £2 



12 



12 



c < fc) - E(t 2 , fe)] - fc) " F(t 2 , fc)]| 



— — sin 2ti V sin 2 ii + c — sin 2i 2 V sin 2 i 2 + c 
6 L 

Collecting the expressions for Jj and using Ji + 3 J2 = 4 J2 — c/ 2 we finally arrive at 

J s = 4J 3 + cli - scl 2 + 4sJ 2 . 



We also need the following integral 



h(h,t 2 ) 



sin * dt 



that evaluates to 



t2 (sin 2 t + c)3/ 2 ' 

h(hM) = -^={[F(t u k)-F(t 2 ,k)] [E(ti,fc)-E(t 2 ,fc)]| 



1 



2(1 + c) 



sin 2* 



sin 2i 2 



V sin 2 *i + c v sin 2 i 2 + c 



3 Menisci 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



In this section we discuss menisci shapes of four types (sphere, catenoid, nodoid and unduloid) between 
two non-touching axisymmeUic solids, sphere and plane. 

First show that for any meniscus type with a fixed sign s and a fixed number n of inflection points 



H^,d x )^H{i>,d 2 ), if d x + d 2 . 



(3.1) 



Indeed, let by way of contradiction, the opposite holds, i.e., there exist two different d\ and d 2 and at least 

one value tp* of filling angle such that H*.(^>*,di) = H ic {^} in d 2 ). However, this contradicts the master 

equation (IA.58I ) which states 

J s (^,^) + nJ 2 (c(^,^))-s(l-cos7rn)/ 2 (7r/2,7r-0 2 ) 
d = — 1 + cos ip* , (3.2) 
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where c(-0*, H+) and I s (ip-k, H+) are given in (12.71 ) and (12. 10b . respectively, and depend explicitly on H+, 
ip*, 6\ and 82, but not on d. The latter means that for these four variables the r.h.s in (13.21 ) reaches its unique 
value which implies the relationship (13- 1 b - 



3.1 Catenoids (c = 0) 

In contrast with the Nod and Und menisci the Cat meniscus exists only for fixed ip values for which H = 0, 
i.e., a = c = 0. We make use of (1A.58I) in the limit c — > in the form (IB.7I ) for s = — 1 where we set the 
l.h.s. to zero to find a relation 

2n — (1 — cos Tm) cos 62 = 0, (3.3) 

from which it follows that the catenoid can be observed only for n = 0, i.e., at the point of the tran- 
sition Undp «-)■ Nod - . Using the relation (1B.4I ) in (1A.58I) for s = —1, n = and approximation 
c = — 4if sin ip sin valid for small H in ( 12.71 ), we find 

2HV = -2HMsinipsinti. 

This leads to 

1 + d — cos ip + sin ip sin(0i + if)) In ^tan ^ 1 ^ — tan -^-^ = , (3.4) 

which can also be derived (see (1A.4I) ) from the YL equation for H = 0. Solutions to the above equation 
exist only for 9\ + if) < it. This condition implies that the logarithmic term in (13.41 ) is negative which leads 
to a stronger condition Q\ + 62 + if) < it. In the special case of ideal plane wetting 62 = this term diverges 
and the equation (13.41 ) has no solutions, so that the catenoids are forbidden. 
Rewrite the equation (13.41 ) in the form 

g{Vl,iP,d) = . . n . + . + In tan — - — = lncot— . (3.5) 

sm(&i + ip) smip sm(t^i + ip) 2 2 

For d = and fixed values of the contact angles g(0i,ip, 0) monotonically grows and tends to asymptote 

at tp = TT — 6\ that implies existence of a single solution of equation (13.4I ). This solution exists when 

the minimal value of g(9i, 0, 0) = mtan(#i/2) is smaller than the r.h.s. of (13.5I) . that leads to condition 

0\ + 62 < 7r. For nonzero d we note that the function g has an additional asymptote at ip = 0. When cZ < 

the monotonic behavior of the function g(6i,ijj,d) does not change, so that we still have only one solution 

to equation (I3.4I ). When d > the function g(0i,ip, d) has a minimum g m j n so that depending on this value 

compared to In cot (#2/2) equation (13.41 ) may have two, one or no solutions. 

Find the maximal value d m of the distance d as a function of ip and 6± for which equation (13.41 ) has a 

single solution. This value is reached when two following conditions are met: dd/dip = and dd/d6i = 0; 
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these conditions lead to 9\ = and we obtain d m = (1 — cos ip m ) / cos ip m , where ip m satisfies the relation 

1 + cos ip m In tan tan — = , 
2 2 



that leads to 



2(1 + d m ) + In 



2 In cot — 



The numerical solution of equation (13.61) for d m (#2) is shown in Figure[2£a). 



(3.6) 




#2 #1 

(a) (b) 

Figure 2: (a) The dependence of the maximal distance d m from the sphere to the plane on the contact angle 
02 for 9\ = 0. (b) The solutions of the equation (13.41) for 62 = 40° for different values of the distance d 
from the sphere to the plane. The dashed curve shows the position of the degenerate Cat menisci. 



There exists a single nonzero value of the filling angle ip satisfying equation (13.41) at d = that leads to 
the Plateau sequence (see Figure [3£a)). Catenoid cannot be observed for d > d m , but there exists a range 
of distances < d < d m for which equation ( 13.41 ) is satisfied for two values of the filling angle as shown 
in Figure He). The corresponding sequence of menisci looks like 



Und r 



Cat — > Nod — > Cat — > Und — > Undj^ — > Undg — > Sph c 



Nod 



The sequence of menisci can pass through a single catenoid but it differs from the Plateau type as shown 
in Figure [Hd). This occurs in the degenerated case when equation ( 13.41 ) has to be complemented by an 
additional requirement: ddi/dif) = 0. Taking derivatives with respect to ip in both sides of equation (13.41) 
and applying the above condition we arrive at a trigonometric equation 



2 sin ip = [ In cot — h In cot — ) sin(#i + 2ip) . 



(3.7) 
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Solution of equation (13.71 ) for fixed 62 determines a curve 0i(ib) on which the degenerated Cat meniscus 
appears. It is shown in Figure |2tb) by dashed curve. 




(c) (d) 

Figure 3: The case (a) presents the typical curve at d = with 9\ = 30°, 62 = 80°, discussed in Q, 
while the case (b) at d = with 9\ = 120°, #2 = 90°, characterized by a requirement 9\ + 62 > vr, was 
not considered in [5]. For positive distance d > one observes (c) with d = 0.4 with 0\ = 0% = 40° 
two catenoids that degenerate for larger distance into a single catenoid (d) with d = 0.53. Different colors 
indicate the menisci of different types: Und^ (red), Undg (cyan), Und^ (green), Nod - and Nod + (black). 
They correspond to colors of types shown in Figure [6] 

3.2 Unduloids (c < 0) 

Before starting to treat the unduloidal solution of equation ( 12.101 ) it is worth to make 

Remark 1 Equation ( 12.71 ) with boundary conditions is associated with nonlinear eigenvalue problem and 
in the case of unduloids its solution has a discrete spectrum and, therefore, is enumerated by two indices. 
The first integer non-negative index n determines the number of the inflection points on the meniscus merid- 
ional profile. When the part of this profile touching the solid sphere is convex the second integer index s 
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takes value ofl, otherwise it equals to — 1. Thus, the unduloid meniscus is denoted as Und^. 

The existence of the Und^ meniscus requires satisfaction of the condition sin 2 t + c > for t G {t2,t{\. 
This condition is rewritten in the form 

sin 2 t > [1 - (1 - 2a) 2 ] sin 2 h . (3.8) 

Denote smt m = min(sin£i, sin £2)- The last condition leads to sin 2 t m > [1 — (1 — 2a) 2 ] sin 2 ii, and 
yields 

a<p-(tp), a>/3 + {4>), where ^ (ip) = ^ (l ± - sin 2 t m / sin 2 tA . (3.9) 

The restrictions (13.91 ) define in the plane {ip, a} a new object B which we called a balloon. It comprises 
an oval and two additional lines a = 1/2 (see detailed description in section [5]). Hereafter the balloon B 
becomes a main tool of our study the dependence a(ip). In regard to H(ip) the balloon undergoes the non 
linear transformation (12.81) . In Figure [4] we present both dependencies H n {i)j) and a n (i/)) with the balloon 
B. 




20 40 60 80 100 120 140 20 40 60 80 100 120 140 



(a) (b) 

Figure 4: The dependencies H n (ip) in (a) and a n (^) in (b) for 6\ = 30°, 62 = 80° and d = 2.3 with two 
disjoined curves. The balloon B is shown by dashed line. Spherical menisci Sph^, Sph^ and Sph^ are 
denoted by black dots. The curve segments correspond to colors of the menisci types shown in Figure [6l 

The relations ( 13.91 ) are independent of distance d and remain valid for all types of inflectional unduloids 
with < a < 1. For all unduloids they define the restrictions on the curvature values that can be obtained 
using the definition (12.8I ). On the other hand, these relations correspond to a single condition on c value 

-sin 2 t m <c<0. (3.10) 

In Figure [5] we present three families of unduloid menisci for fixed filling angle ifi (Figure Sa)), volume V 
(Figure |3b)) and surface area S (Figure [3c)) of pendular ring. 
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(a) (b) (c) 

Figure 5: The menisci Und^ (red), Und^ (green), and Und^ (blue) for 9\ = 30°, 62 = 80° and d = 2.2 
computed in three different setups: (a) for %jj = 60° they have different volumes V and surface areas 
S; (b) for fixed V = 3.6 there are 4 menisci with Vi = 34° and ip 2 = 47.2° - Und^, V>3 = 61.7° - 
Und^, ^4 = 71.3° - Und^~ and different surface areas S; (c) for fixed S = 11.3 there are 4 profiles with 
tpx = 30.3° and ip 2 = 47.1° - Und^, ^3 = 56.7° - Und+ ip 4 = 61.95° - Und+ and different volumes V. 

3.3 Nodoids (c > 0) 

The Nod meniscus has positive c irrespectively to its concave or convex version. Therefore the above 
analysis of a fails and it is replaced by another one, less strong but still universal. 
Show that for d > the Nod s menisci satisfy the following constraints, 

H(ip, 0) < H(ip, d)<0 (Nod"), < H(ip, d) < H(tp, 0) (Nod+) , (3.11) 

using (13.11 ) and several additional observations listed below. 

First, the curvature H for any finite d is always negative for the Nod - meniscus and always positive 
for the Nod + meniscus. Next, the Nod - meniscus disappears at finite d when two catenoids annihilate. 
Finally, when %jj — > it and d — > 00 the curvature of the Nod + meniscus tends to zero. Combining these 
facts with (13.11 ) we arrive at constraints (13.1 lb . 

The main statement stemming from (13.111 ) is that the two different curves H(tp,di) and H(tp,d2), 
d\ ^ c?2, do not intersect in the entire angular %jj range of the Nod - and Nod + menisci existence. Although 
this statement is of high (topological) importance, it does not provide the quantitative estimates. Therefore 
for the Nod - meniscus we give one more estimate for the curvature. 

Start with relationship between the curvature H and the surface area S for Nod - meniscus. For this 
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purpose combine formulas (l2.10[|2?Tll,l2.13ll2.14l|2.19l ) and obtain 

2H 2 



-S = K-=2(h-h)-h = WV-h. (3.12) 

TT 

Keeping in mind the positiveness of S we arrive at the bound AH^f > I3. Using the definition (12.161 ) we 
obtain ^3 > —b^fc, where 5 = tt — 0\ — 62 — ip, 5 > 0. Thus, substituting this relation in the above 
inequality we find 

H>-^V~c. (3-13) 

This bound does not contradict inequalities ( 13.111 ) for the Nod - meniscus since its curvature is negative. 
Keeping in mind this fact and substituting (12.71 ) into (13.131 ) we arrive at 



sin ip sin(#i + ip) 
sin 2 ip - 4^/ 2 5~ 2 



H > . I , \,t, r _o • (3-14) 



a > r0 . 2 y . (3.15 

t z 1/1 _ A \T;^ 



This inequality is equivalent to 

5 2 sin 2 ip 
5 2 sin 2 ip - 4^ 2 

The last inequalities have one important consequence. Since H < then the necessary condition for the 
Nod~ meniscus existence is 

5sinip <2<S/ - — 0l - ° 2 — -sinip + cosip <l + d. (3.16) 

For convex nodoid Nod + meniscus it is possible to find explicit expression for the upper bound H(ip, 0). 
First, note that as it is demonstrated in Appendix [C] meniscus curvature has a local minimum, so that the 
upper bound is given by the largest of two curvature values - H{<Pq) at the sphere Sph,^ and H(tt) at 
ip = n. These values are 

W0) sinc/>+ ' K ' 2 + d 

Show that H(<Pq) > H(tt). Indeed, this condition implies 

(1 + d) sm{9\ + 4>q) + cos 62 sin^Q - > sin^ — sin(0i + 0q ). 

Using in the above relation the l.h.s. of formula dB.9b we obtain 

sin 6\ > sin^Q — sin(#i + <^), 

leading to 



cos — + cos U$ + —J = cos cos -y > 0. (3.17) 
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Recalling that Q\ + (f>^ < it we see that (13.171) is always valid. Thus, the upper bound for the Nod + 
meniscus is given by H(cf>Q , 0). To find this value explicitly we use dB.9b with n = and d = to produce 

sin(6»i + 4>q) + cos 2 sin $ - sin 6\ = =>• cos ( Q\ + J + cos 2 cos -y = 0- 



The last relation implies 



COS 01 + COS 6» 2 

tan-^- = — . (3.18) 

2 sin 0i 



On the other hand, we have 

0) = sin (^ ±#l) = cos 01 + sin ^ cot ^+ . 
sin 0q 

Using (13.181 ) in last expression we find for the upper bound of Nod + meniscus 

H(ti,0) = - f 2 * 2 (3.19) 

2 (cos 01 + COS 02 J 

3.4 Spheres (c = 0) 

The spheres can be considered as a limiting case of the unduloid menisci and, therefore, also labeled by 
two indices Sph^. Consider a function 

f n (tP) = (1 + d) sin(0i + V) - (n - -5) sin ^ - sin0 1; (3.20) 

which root ip = <ft+ satisfying the equation fnifin) = provides the value of the filling angle at which 
sphere Sph+ is observed (see (1B.9I) ). In (13.201 ) one has S = for odd n and 6 = cos 02 for even n. It is easy 
to check that 

f n +l(i>) ~ fn{4>) = -[1 + cos 7m cos 2 ] sin^ < 0, 

which means that the curve / n +i lies below the curve f n . 

The value of the function f n at %jj = reads / n (0) = <isin0i which is non- negative, while f n (ir) = 
— (2 + d) sin 0i is always negative. Noting that the function /„ is a periodic one with the period 2tt we 
find that for positive d in the interval {0, tt} the function /„ has a single root and f n {4>n) < 0. As 
f n+ i(4>+) < f n {4>t) = and /n+i(^«+i) < 0we immediately find that < </>+. This means that for 
positive d the value tp n of the filling angle at which sphere Sph+ is observed decreases with increase of n. 
For very large n the value of </3+ tends to zero. Using (13.201 ) we find in linear approximation 

■ a+ a+ dsaxOi dsin0i 

sm^ = 0n= r , n»l. (3.21) 

n — o — cos 0i n 

For d < we have / n (0) < and the first derivative reads f' n (0) = cos0i — (n — 5). For n > this 
derivative is negative, so that the function /„ has no roots. When n = we have /q(0) = cos 0i + cos 02- 
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This expression is negative for 81 + 82 > n, so that no spheres exist for d = when the last condition holds. 
Finally, in case of the wetting sphere 8\ = we have f n (0) = f n {^) = and no spheres are allowed to 
exist. 

Show that for fixed n the value (j>n{d) of the filling angle of sphere Sph+ increases with growing d. 
Indeed, from (13.201) one finds that the difference 

f n ( ( f>t(d),d 1 )-f n (^(d),d) = (d 1 -d) S m(8 l + ti(d)), 

is positive for d\ > d. As f' n {4>n) < it immediately follows that > <f)n{d). 

Another type of sphere Sph^ (n > 0) discussed in IB .21 can be considered as a limiting shape of the 
Und~ meniscus at ip — > 0, its curvature is given by (1B.8I) . 

In Table 1 the characteristic signs of H, c are given for different types of menisci. 



Table 1. 





Nod" 


Cat 


Und s n 


Cyl 


S P h s n 


Nod+ 


H 







+ 


+ 


+ 


+ 


c 


+ 












+ 



Note that the meniscus Und„ comes in several different types depending on the number n of inflection 
points and the curvature of the segment of meridional profile is touching the sphere: s = — 1 for concave 
profile and s = +1 for convex one. 

4 Unduloid Menisci Transitions 

Transitions between unduloids and nodoids of different types can be easily classified. Namely, there exist 
transitions between the concave Und^ (convex Und^) unduloid to concave Nod - (convex Nod + ) nodoid 
through catenoid Cat (sphere Sph^). The transitions between different types of unduloids are numerous and 
schematically presented in Figure [6] As shown in IA.7l and lA.8l all inflection points for a given inflectional 
unduloid Und„ have the same abscissa value = smtl/(2H) = \/^c/(2H). An addition (removal) of 
an inflection point may take place only as the result of the inflection point separation from (merging with) 
the sphere or the plane. When the inflection point is on the sphere we have c + sin 2 t\ = to obtain using 
(TAJl) 

2Hsmip = sinii. 
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n A 



Sph 3 



Und; 



Undt 



SphJ 



Sph 2 " 



Und; 



Undi 



Sph+ 



sphr 



Und7 



Undl 



Spht 



Nod" 



Nod + Cat 



Undn 



Cyl 



UndJ 



Nod + 



Sph+ 



1/2 



Figure 6: Possible transitions between menisci of different types shown as a function of the parameter a. 
The left vertical dashed line represents the transitions at a = (3~ between Und~ unduloids. Similarly, the 
right vertical dashed line corresponds to the transitions at a = f3 + between Und+ unduloids. Finally, the 
transitions between unduloids of opposite signs takes place at a = 1/2 and is represented by the central 
vertical dashed line. 



Recalling definition (12.81 ) we find a = 1/2. Noting that menisci Und^ (Und+) exist at a < 1/2 (a > 1/2) 
we find that transitions of the type Und„ <-> Und^ occur at a = 1/2 when the inflection point is on the 
sphere. For the inflection point on the plane we have c + sin 2 t2 = to obtain with (IA. 17b 



2H sinip = sinti + s\/ sin 2 ti — sin 2 12 , 

and we find using ( 13.91 ) that in this case the transitions of the type Und„ «-)■ Und;^ occur at a = (3 s . 

Menisci of the Und+ type exist for 1/2 < a < 1. With growth of a the abscissa of the leftmost point of 
the meniscus meridional section tends to zero and reaches it at a = 1. At this moment the profile is made 
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of several segments of a circle and the meniscus touches the axis of rotation x = 0. The spherical menisci 
Sph+ata = lare described in lB.ll 

The above considerations lead to the following rules of transition between unduloids: 

Und s n <— ► Und s n±1 at a = /3 s , 
Und s n < — > Und"^! at a = 1/2. 

Thus we show that the balloon introduced above can be viewed as a set of transition points between undu- 
loids. 

Show that for fixed values of Q\ , 8 2 and d the curves corresponding to different unduloid types never 
intersect (except for the transition points discussed above that arises in case when the unduloid orders differ 
by unity). As the curves corresponding to unduloids of opposite signs cannot intersect in {tp, a} plane, we 
have to consider only unduloids of the same sign. 

Consider first same sign unduloids of the orders that differ by an even number, for example, Und| k and 
Und| k ,, or Und| k+1 and Und| k , +1 , where k 7^ k'. It immediately follows from (1A.46I ) and (IA.52b that the 
difference between the curvatures of these menisci is {k — k')I 2 /^> 7^ 0. 

If the order difference is odd and larger than two we have for unduloids Und| k and Und| k , +1 with 
k' 7^ k, the curvature difference reads [(k' — k)I 2 + ^2(^*5 ^2)]/^- Consider the integrals I 2 (t^,t 2 ) which 
are finite and real. Note that i+ < ti < t~ , which leads to 

|/2(02)| <h(t-,t+) = / 2 - 

The last relation implies that the curvature differences mentioned above are nonzero that finishes the proof. 

From (IB.7b it follows that the curvature of spherical menisci Sph+ grows monotonically with order 
increase without any restriction to the order value. This means also that the corresponding unduloids Und+ 
can be observed without any restrictions to the order n. 

The asymptotic behavior of the curvature at small filling angles ip <C 1 is discussed in Appendix IB.ll 
For d > at zero filling angle ift = only concave unduloids Und~ (n > 0) exist. From (IB. 8b it follows 
that for n = the curvature H (0) = leading to contradiction as the curvature can turn to zero only for 
catenoid. It is important to underline that there are no other restrictions to existence of Und~ (n > 0) 
unduloids at zero filling angle. 

Using the formula (1A.58I ) for unduloid curvature we have 

2^* = 2i(ti,t 2 ) + sl 2 (ti,t 2 ) +nl 2 - s(l - cosvrn)/ 2 (7r/2,t 2 ), a = sgn(2a n - 1) , (4.1) 
where the sign s = ±1 corresponds to two menisci which differ by the sign of the meridional curvature 
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at the meniscus-sphere contact point. Equation (14.11 ) defines the function a n (ip) in two different regions, 
a n > 1/2 and a n < 1/2, while a n (ip) is a smooth at a n = 1/2. 
Rewrite (14.11 ) as follows, <3?* (a ra , ?/>) = 0, where 



$ s n (a n ,ip) = -2a n *sinti + Ii(h,t 2 ) + sl 2 (h,t 2 ) + n/ 2 - s(l - cos irn)I 2 (ir/2, t 2 ) 
One can define the derivative o/ n {ip) having a unique value determined from the equation 



sin-^, (4.2) 



A s n aM + B s n = 0, A s n 



DS 

da' n 



5-0 ' 



(4.3) 



where the both functions ^4* (a n , ^) and B^(a n , ifi) do not vanish simultaneously. 
Direct computation gives the following expressions for A s n {a n , ip) and B^(a n , ip) 



A s n {a n ,ip) 



B s n {a n ^) = 2a n V 



-2^sinti - 

2(2a n — 1) sin 2 1\ sin ip 

sin 6>i 4a 2 - 2a n - 1 . . 2 

sin ti sin w + sin w — 

y 2<t„ — 1 



sl/i(ti,t 2 ) — 2nl 2 (c) — s(l — cos 7rn)/4(7r/2, i 2 ) , 
i 



(4.4) 



sini/; 2a n — 1 

2a n (a n — 1 ) sin 2t i sin ^ 



5/4(^1,^2) - 2nl 2 (c) - s(l - cos 7m)I i (TT/2,t 2 ) , (4.5) 



where the integral h(ti,t 2 ) is computed in (12.241 ) and the derivative I' 2 (c) is given by (ID. 13b . 
4.1 Transitions Und~ <h- Und+ +S 

Consider first the transitions on the line a = 1/2 between unduloids of opposite signs. This transition takes 
place when the point t\ = separates from the sphere, where t+ = arcsin c and t~ = it — arcsin \/— c. 
As H = sin ii/ (2 sin ip) we use d4.il) to obtain 



^i(**,*2) " I* 2 {tlM) + nh + (1 - cosvrn)/ 2 (^/2,t 2 ) sin(^ - 0! 



(4.6) 



where we introduce a special case of the integral I 2 



sin 2 i dt 



t 2 v sin 2 i — sin 2 ti 



(4.7) 



The general expressions for the abscissa of the meniscus meridional profile contain the term \J sin 2 t + c, 



which at the transition point transforms into y sin 2 i — sin 2 t^r. It leads to a condition t 2 < producing 

t 2 <n-ti — > ip<6 2 -8 ll s = +1, 
h<h — > tp >tt-6 2 -9i, s = -1, 

implying that the transition takes place to the left (right) of the balloon for s = 1 (s = —1). 
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Show that the transition considered in this subsection is smooth, i.e., unduloids Und~ and Und^ +S 
meet smoothly at a = 1/2. It means that the value of the derivative a'(tb) computed on both sides of the 
transition point is the same. As at the transition point we have c + sin 2 ii = the integral l4,(ti,t 2 ) 
in (14.414.51 ) diverges. We show below that nevertheless both A and B have finite value at the transi- 
tion point. In the vicinity of a = 1/2 we introduce a = 1/2 + se and find c = (4e 2 — l)sin 2 ti and 



yc + sin 2 t\ = 2esin£i. Making use of (12.241) decompose the integral h(ti,t 2 ) into diverging Iu{ti,t 2 ) 
and non-diverging I& c {t\, t 2 ) parts, 

h{h,t 2 ) = hc(ti,t 2 ) + hd(ti,t 2 ) , hd(ti,t 2 ) = j- n (4.8) 

s{2a — 1) cos ti 



Using this relation in (14.41 ) we find 

2 sin^ £i sin ib 

A s = -2^> sin ti = 2 tan ti [cos 6 1 - (1 + d) cos h] . (4.9) 

costi 

Substitute (14.81) into (14.51 ) and note that two diverging terms cancel each other in vicinity of a = 1/2, 

4a 2 -2a — 1 , . , . , 

— sinti sin^y — 2sa(a — 1) sin 2t\ sinwi^fii, t 2 ) 

la — 1 

1 1 sin 2ti sin^ 

sin t\ sin tp = 0. 



2a -1 r 2(2a-l) cosh 

The remaining terms read 



B^ = ^ p + sin 2 tp + sin t\ cos t\ sin tp 

smip 



sl4c(h,t 2 ) — 2n/ 2 (c) — s(l — cos 7rn)/4(7r/2, t 2 ) • 

Show that this expression is conserved at the transition Und~ -R- Und^. It is sufficient to show that it is 
valid for the expression in the square brackets in (14.1 II ) that leads to the relation 

hc{t\M) ~ h(w/2,h) = sii' 2 (c), 

where s\ = 1 (s± = — 1) corresponds to the transition to the left (right) of the balloon with t\ = t* 1 . The 
last relation can be written as follows: 

hM\7r/2) = Sl i 2 (c), (4.10) 

It is easy to see that 

hc(tt,7T/2)=h c (7T/2,TT-tt) = il 4c (^,7T-^), 

and using ( ID. 19b we establish the validity of (14.10b . Thus we find at a = 1/2 

B s = — i + sin 2 tp — sinti cos t\ sin tp[(2n — ssi)/ 2 (c) — s cos irn I An 12, t 2 )\. (4.11) 
smtp 

Selecting here s = — 1 we arrive at the final expression for B 

B~ = ^ Sm ^ 1 + sin 2 tp - sinti costi smip[(2n + si)/ 2 (c) + cos7rn IJit/2, t 2 )], (4.12) 
sin tp 

where si = 1 (si = — 1) corresponds to the transition Und~ <-> Und+ + to the left (right) of the balloon. 
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4.2 Transitions Und* <-> Und^ +1 



This transition for odd n takes place at a = (5 s = (1 + syl + c/ sin 2 ti)/2 when the point i 2 = 
separates from the plane and we obtain 



2Htf! = h(t u t-') + a/ 2 (ti.O + (" + 1)^2. 



(4.13) 



As c = — sin *2 we find 



sin 2 t 



i2 



sin 2 1 



and arrive at 



t 2 \J sin 2 t — sin 2 i 2 Jti \J sin 2 i — sin 2 T 2 

sin^ 



-ij(t2,*l), 



(4.14) 



/i(ti,t-')-s^(Cti) + (" + lK2 



2/? s sinii 

For even n this transition is observed with separation of the point t 2 = t% from the plane for which we 
obtain 

r _ / . . , . . ± l sin ib 

L v * y 2V * y J 2/3 s sinti 

Discuss one more question: how smooth are transitions Und+ <-> Und r [ +1 at upper (a = /3 + ) and Und~ -H- 

Und~ +1 at lower (a = (3~) arcs of balloon, respectively. To this end, consider equation (14.31 ) in vicinity 

of the transition point {ip*, f3 s (ib*)) belonging to the balloon and estimate the leading terms in (14.41) and 

(14.51) when c — > — sin 2 t 2 . The only divergent terms are integrals I±{t\, t 2 ) and 1^/2^2) which by (12.241 ) 

behave as follows, 

r+sin 2 fri — » n sin 2/o 

(4.15) 



T / \ T / /„ \ c+sin t2 — > 



2(1 + c) A/sin 2 t 2 + c 

Substituting (14.151 ) into ( 14.41 ) and ( 14.51) we arrive in the limiting case c + sin 2 t 2 — >• to explicit formula 



(4.16) 



2/3 s - 1 

which is finite for f3 s ^ 1/2 and independent on order n, i.e., both curves Und^ and Und^ +1 meet smoothly 
at the balloon. 



4.3 Transitions Und 1 -h- Und^ , Und 2n Und^ and Und 2n <-»■ Und2 n+2 

In previous sections we have discussed in details conditions and rules of regular transitions between un- 
duloids that accompanied by addition or removal of one inflection point. It is also possible to observe 
degenerate transitions when the number of inflection points changes by two or does not change at all. Find 
the conditions for such degenerate transitions. 
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An addition of two inflection points takes place only when one of these points separates from the solid 
sphere and the other one from the plane. The first event corresponds to a condition c + sin 2 t% = 0, 
while the second one requires c + sin 2 £2 = 0. These relations imply sinii = sin £2 leading to ip\ = 
min{#2 — 0\,tt — 6\ — 62] and ip\ = max{# 2 — #i,7T — #i — 62} ■ These critical values define the extremal 
points of the balloon where it transforms into the segments of the line a = 1/2. When ip\ is negative only 
a part of the balloon is observed. When both values ip* are negative the balloon does not exists. Finally, for 
62 = vr/2 the balloon reduces to a point at ip* = tt/2 — B\. 

Consider the transition Und^ Und^ k+2 at ip = ipf. Using (IA.46I ) and noting that t\ = 62 = t~ we 

find 

= h(9 2 ,ir- 2 ) - h(92,TT- B 2 ) + 2kl 2 = -2 cos 2 + (2k + 1)J 2 , 

simp 

from which we obtain the distance d n where this transition occurs 



dn = -l + cos(g 2 - fli) + Sm(g2 - 9l) (n+l)/ 2 -2cos( 

sin 62 



(4.17) 



The case k = should formally correspond to the transition Undg <-> Und^ - It can be checked that Undg 
unduloid can exist only below of the balloon, and not to the left or right of it, so that the above transition is 
forbidden. The relation (14.171 ) at k = leads to a special case Und^ «-)■ Undg Und^ when the segment 
corresponding to Undg reduced to a point. The sequence of distances given by (14.171 ) is a periodic one with 
the period equal to 

, 2/ 2 sin(# 2 -fli) , . fa a . w {, 1 
d n+2 -d n = — n = 4sin(0 2 - Bi)E 1 



sm U2 V sm V2 , 

It can be shown that the transition Und^^ Und^ k+3 is forbidden at ip = ip* . As the number of inflection 
points is odd both new points to be added should correspond to either t+ or t~. It means that t\ = i 2 which 
contradicts (for 62 7^ vr/2) to the ip\ value. 

The second critical point ip = ip\ initiates the segment of the line a = 1/2 to the right of the balloon. 
Only the transition Und~ «-)■ Und,}^ is allowed on this line when the inflection point merges the sphere. 
It means that at ip = ip*, the inflection point with t = t% merges the sphere while the point with t = i 2 
separates from the plane. As the result total number of inflection points remains constant. For the Und^ 
unduloid we find t\ = i 2 that corresponds to ip = ip*,. Thus, the transition Und^ 44> Und^ is allowed, 
while Und2 k+1 44> Und^ k+1 is forbidden again (for 2 7^ 7r/2). In this case we have 

sin 82 



sin ip 

from which we obtain 



* = 2kl 2 , 



-l-cos^ + ^+n ^ + ^ J.. (4.18) 

sm O2 
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The sequence of distances given by (14.181 ) is also a periodic one with the period equal to 

dn+2 - d n = 2 A^i±M = ism( o 2 + e x) E ( i - -1 



sin 9 2 V sm #2 

5 Topology of Unduloid Menisci Transitions 

This section is mostly topological and deals with qualitative behavior of different branches of function 
a(ip) enumerated by the number n of inflection points at corresponding menisci. We list the most general 
properties of curves a n (ip) in the plane {ip, a} and study ramification of these curves around the balloon. 
This global geometrical representation allows to classify possible trajectories and their intersection points. 

5.1 Balloons ^{ip) 

Consider the rectangle A := {— 9\ < ip < tt — 9i, < a < 1} in the plane {ip, a} and define a balloon B 

M = B l \JB r \JB d \JB u , where (5.1) 

Bl = j-01 <1p < @min , OL = i| , B d = {Qmin < Ip < Qmax , « = /3~ (ip)} , 

B r = l@max < i> < K ~ 9 X , OL = , B u = {Qmin < 4> < Qmax , « = P + {ip)} , 

and Qmin = min{02 - 9\, ir - 9 2 - #i}, & max = max {#2 - 9\,ix — 9 2 — 9{\, < 9\,9 2 < it. Two 
functions fi + (ip) and f}~ (ip) give the upper and lower parts of convex symmetric oval, 



Subscripts /, r, d and u stand for the Ze/?-, right-, down- and upward directions on B. Denote by B the 
balloon B with its open interior 25, 

F = BU<B, ® := {9 min < V < 6^ , <«</3 + (^)} • (5.2) 

In special case 9 2 = vr/2 we have B^ = {— #i < ip < 7r — #i , a = 1/2} while the part Bd U 5 U of 
balloon is reduced into a point O = {ip = n/2 — 9%, a = 1/2} such that G B^/2- 

5.2 Trajectories a n (-0) 

Below we give a list of rules for topological behavior of a n (ip) in the presence of B. 
1. a n (ip) is completely defined by three parameters: < 9\,9 2 < ir and d > —2. 
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2. a n (ifi) is a real function representable in the {ip,a} plane by a nonorientable trajectory Y without 
self-intersections. All trajectories are continuous smooth curves and located in domain A' \ 53, where 

A' = {0 < ip < 7r - 0i, < a < 1}. 

3. The final points of F are associated with spheres Sph^ and Sph+. Equip T with indices according to 
the final spheres designation in such a way that a left lower index does not exceed a right lower, i.e., 
ni r n 2 := {Sph~ — > Sph+}, ni < n 2 . The following coincidence property holds: ~T+ = +T~. 

4. In the {tp, a} plane the spheres satisfy 

Sph- G {if) = 0, a n = 0} , Sph+ £ {0 < ip < tt, a n = 1}, 01 > , (5.3) 
Sph- G {V = 0, < a n < 1} , Sph+ G = vr, < a n < 1}, 0i = . 

5. Angular coordinates %jj = of spheres Sph+ are arranged leftward in ascending order (see IB .3b . 

< . . . < 0+ < < . . . < (f)^ < 4>^ < tt - 6»i . (5.4) 

6. There exist generic trajectories of four topological types, 

n^n+ii n^n+l' n— l^n > n-l-^n i n> 1 . (5.5) 

7. Different parts of trajectories are labeled by different sub- and superscripts Und^ and Und~ 2 where 
the upper index is equal to s = sgn(2a n — 1). 

8. In vicinity of the point (tp = 0, a n = 0) the sheaf of trajectories „T~ +1 , „T+ +1 and „_iT~ is 
build in such a way that slopes £~ of the Und~ parts are arranged clockwise in descending order (see 
section IB .21 ). 

| > . . . > c > e r ;_i > . . . > > > ° • ( 5 - 6 ) 

9. For given 0i, 62 there exists a unique d = such that there appears one of five types of intersection 
(saddle) points: 

a~ — ~r+ 6?> + r~ a+ — ~r+ a ~~ r+ 

"n — n^n w n—l'~n+l ) "n — n 1 n w n-l 1 n+1 ' 

B n = n^n ® jt_iT^ +1 , B+ = ra _ir+_ 1 ® n-2^n > = „ (g> ^_ir+ +1 , (5.7) 

where operation Ti ® T2 denotes intersection of two trajectories Ti and T2. The indices of a saddle 
point correspond to unduloid Und^ observed at this point. The saddle points B° of mixed type are 
located on a line a = 1/2 (s = 0). 
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(a) (b) 
Figure 7: Plots a n (V>) for (a) 6>i = 10°, 6 2 = 70°, d = 6.45 and (b) B x = 10°, 6 2 = 70°, d = 6.7. 

10. Every trajectory of the types „T~ +1 , „T^ +1 and ^_ x r~ is necessary smooth at a = 1/2 (see 
section 14. II ). 

11. Every trajectory of the types „T~, x and n^n+l * s necessary tangent to Bd U i3 u (see section |4~2~1) . 

12. The changes of indices in Und*. occur at balloon B in accordance with Table 2 (see Figure 13). 

(a) The change of the upper indices (s «-> — s) and lower indices (ni <-> n 2 ) occurs at the point 

(b) The change of the lower index (m n 2 ) only occurs at the point T € £>d U B u . 
In the nondegenerate case one has |m — ra 2 | = 1. 



Table 2 





n 1 n+1 


n n+1 


+ r- 

n— 1 n 


+ p+ 
n-1 1 n 




— 7""+ — 7+ 
n J «+l' n+l- t n+2 


- T + 




- 7-+ 

n— 1 n 








+ T - 
n— 1 n 


+ 7 1 — +T~ 
n— 1 n ' n ^n+l 


fid 


n-'n+l 




~~ T~ 
n— 1 n 


n ^n+l 


fi« 


+ T + 

n+l ± n+2 


+ T + 
n+l J n+2 




+ 7^+ 
n—ln 



In Table 2 a symbol denotes a point belonging to two parts Und^ and Und,[ 2 of trajectory. Empty 

boxes mean that corresponding transitions do not exist. The trajectories can be tangent to balloon Bd U B u 
at its left and right points (see Figured]), 

r T i + » 2 -r+ +2 e ^n{H d UH u }, and 2 „r+ e fi r n{fi d ui3 u } . (5.8) 
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6 Saddle Points 

For given values of the contact angles 61,62 a change in the distance d value leads to changes of the 
trajectories shape in {ip, a} plane. Sometimes such transitions are accompanied by drastic changes of the 
trajectories' topology characterized by an appearance of the saddle points. The saddle point can be defined 
as a point that belongs to two trajectories simultaneously. 

It is instructive to find the coordinates of the saddle point {ip^, a^} as well as the distance d c n at which 
the saddle point is observed. Below we describe a procedure for such computation for each type of the 
saddle points belonging to a single meniscus type Und„. 

6.1 Saddle Points of Simple Type 

First note that the saddle point may be observed at the intersection of two segments of the curve a n (ijj) 
determined by the sign s and order n of unduloid meniscus Und^ that completely defined by the relation 
(14.21) . At every point of these segments (except the saddle point) one can define the derivative a' n (ifj) having 
a unique value determined from the equation ( 14.31 ) where both A s n and do not vanish simultaneously. At 
the saddle point the derivative a4(V0 is not unique and in this case A s n = B^ = 0. Thus the saddle point 
{^,a^}at(i = (ijjis determined from the condition A s n = B^ = together with (I4.2I ). 
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(a) (b) 

Figure 9: Plots a n (ip) for 6 1 = 30°, 9 2 = 90° at (a) d = 7.66209 and (b) d = 5.82697. In (b) a vicinity of 
the saddle point is shown. 



Using the condition A s n = we find the expression in the square brackets in (14.41 ) and substituting it 
into (14.51) we obtain at the saddle point 

„. sin#i 4a n (l — a n ) T 4a 2 — 2a n — 1 , 9 , 

5* = 2a„tf i p ^tfcosii 2 sintisinV + sin 2 ^. (6.1) 

sin^y 2a n — 1 2a n — 1 

Use (16.11 ) in the condition = to express \f at the saddle point 

sin 2 if) [(4a 2 — 2a n — l) sinii — (2a n — 1) sin V>] 



and eliminate 



2a n [(2a n — 1) sin#i — 2(1 — a n ) costi sin^] ' 



d c n = ¥£ + cosV£-l, 



(6.2) 



(6.3) 



from the saddle point conditions. Thus we arrive at the final equations determining the saddle point position 
in the {ip, a} plane: 



2a n \I>^sinii - Ii(ti,i 2 ) + siz(*l>*2) + n h - s(l - cos 7rn)/ 2 (7r/2, i 2 ) 



sinV = 0, (6.4) 



^ + (2a n — 1) sinii simp sl^iti, £ 2 ) — 2ral 2 (e) — s(l — cos im)!^ /2, i 2 ) 



0. (6.5) 



Solving the equations (16.416.51 ) we find the saddle point. An example of trajectories in a vicinity of such a 
point is shown in Figure [TOf a). 



6.2 Saddle Points of Mixed Types 

The saddle points considered above belong to a single type of unduloid meniscus and are characterized, in 
particular, by the sign of a n — 1/2. When at the saddle point a n = 1/2 this point belongs to two different 
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(a) 



(b) 



Figure 10: The saddle points at intersection of ~T+ and n-i^n+i trajectories for n = 1. (a) The trajectories 
in the vicinity of saddle point for 6\ = 30°, 62 = 60° and d = 1.05187. (b) The trajectories in the 
vicinity of B° mixed type saddle point observed for 6\ = 30°, 62 = 80° and d = 2.2465. Colors and 
numbers correspond to indices of respective unduloids. 

types of menisci of different orders and signs. These points do not have unduloid sign characteristics and we 
designate them by a smaller (of two) order only. Consider first such saddle points located to the left of the 
balloon (Q m i n , @max > 0). This transition described in Section l4~T1 happens when t\ = t+ = arcsin v / — c. 
Equation ( 16.51 ) cannot be used directly at a n = 1/2 as integral I4 diverges at t\ = tf. In the limit a — > 1/2 
we have by (12.24b . 



1 



and obtain from (16.5 



|1 - 2a| % /TTc 

= sin^tant+. (6.6) 
Note that this value does not depend on both order n and contact angle 62. Substitution of (16.61 ) into (16.41 ) 



produces a condition on c value 



h(tt,t 2 ) - h{tt,t 2 ) + nl 2 + (1 - cos7rn)/ 2 (7r/2,i 2 ) + 



0, 



(6.7) 



where order n corresponds to the meniscus with a n < 1/2. For given values of order n and contact angle 
62 we find c value verifying the last condition that leads to determination of ip* = arcsin \f^c — 6%. Using 
it in (16.61 ) and (16.31) we arrive at the distance value for which the mixed saddle point is observed. An 
example of such a point is shown in Figure [TOl b). 

Consider a mixed type saddle points located to the right of the balloon. This transition described in 
(14.11 ) happens for s = — 1 when t\ = t~ = it — arcsin y/^c. Equation (16.51 ) leads to 



^ = sin ijj tan . 
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Using it in (16.41 ) we obtain a condition on c value 

h(tZ,t 2 ) ~ h(t-,t 2 )+ni 2 + (1 - cos7rn)/ 2 (7r/2,t 2 ) - -£= = 0, (6.8) 

Vl + c 

where positive order n corresponds to the meniscus with a n < 1/2. Show that the equation (16.81 ) does not 
have solutions for n > 0. Using the relation 

h{t-,t 2 ) = y + h(n/2,t 2 ), 

rewrite the left hand side of (16.81 ) as 

, + cos £2 + (n — 1/2)22 — cos7rn/2(7r/2,t2), 
yl + c 

where sum s\ of the first two terms is always positive. For odd n = 2k — 1 we have s\ + (2k — 1/2)I 2 + 
l2(7r/2, i 2 ) > 0. For even n = 2k we have s 1 + (2fc - 1/2)J 2 - I 2 (ir/2,t 2 ) > 0, as I 2 (n/2,t 2 ) < 
I 2 (ir/2,t^ : ) = I 2 /2. Thus the saddle points of the mixed type cannot be observed to the right of the 
balloon. 

6.3 Saddle Points Sequences 

The computation of the saddle point position for fixed values of the contact angles 6\ , 9 2 and increasing n 
shows that for large n a sequence {ipni a n} accumulates in a small vicinity of a point {ip%, a£} belonging 
to the balloon (not reaching it), i.e., a% = /3 S (^) as shown in Figure [TT1 

It is instructive to determine the position of the accumulation point {tp^, a^}. First note that in ( 16.41 ) for 
n»l the dependence of ^ on n is determined by a relation 

2a^^sinti = nl 2 , 

implying that both ^ c n and d c n grow linearly in n. As the integral I' 2 (c) is always negative for c < then 
the leading term for n ^> 1 in ( 16.51 ) is contributed by the integrals 1^. From ( 12.241 ) it follows that the proper 
divergence is provided by the term sin2t2/[2(l + c) y sin 2 t 2 + c] when c — > — sin 2 t 2 . Substitution of 
(3 s (ip) in ( 12.91 ) shows that c + sin 2 9 2 = holds on the balloon. Thus the divergence of ^ used in (16.21) 
leads to the condition 

[2at - 1) sin 5i = 2(1 -Osin^cos^ + V, c ), (6.9) 

where 

c if, , L ^2 
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Figure 11: (a) The saddle points (green) and A^ k for 9\ = 30°, 62 = 60° and k = 2, 3, 4, 7, 10 approach 
the accumulation point (black) on the balloon with growth of the index 2k marking the saddle points. The 
inset shows saddle points (cyan) for k = 1, 2, 3, 4, 7, 10 approaching the accumulation point, (b) The 
saddle points (blue) h^k-i ^ or ^1 = ^0°, 02 = 120° and k = 2,3,4,5,8, 11 approach the accumulation 
point (black), (c) For 6>i = 30°, 2 = 88° the saddle points B~ for n = 1,3 (red), n = 2,4,6,10,16 
(cyan) and B+ (green) for n = 6, 8, 10, 16 approach the corresponding accumulation points (black), (d) 
For 0i = 30°, 2 = 90° the saddle points B" for n = 1, 3, 5, 7 (red) and n = 2,4, 6, 8, 10, 14, 20 (cyan) 
approach the accumulation point (black). 

It follows from ( 16.91) that the accumulation point with a% < 1/2 can be observed for > tt/2 — Q\, while 
at the point with 01% > 1/2 we have < 7r/2 — 0\. The relation ( 16.91 ) also implies 

2(1 - a%) sin(0i + tf) cos = sin X , 2a% sin(0j + cos ^ = sin(0i + 2$;). 



Multiplying these relations and using the definition ( 12.91 ) we find the sign independent condition on the 
accumulation point 



sin 0i sin(0i + 2^) = sin 2 2 cos 2 ip c 



(6.10) 
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It can be shown that the equation (16.101 ) has two solutions corresponding to the accumulation points 
belonging to (3 s (tp). Computing ip° from (16.91) or (16.101) we find a growth rate of d c n for large n which is 
given by 



sin tpl 



-h (- sin 2 e 2 ) = - . S Z 2 ^\ IC J 2 (- sin 2 e 2 ) . 



2a% sin (9 1 + ifc 



2 sin (0i + 2V£ 



7 Touching and Intersecting Bodies 

When the distance d between the solids is non-positive it leads to strong simplification of the topological 
structure of the solutions of (I2.ll ). First, only a single branch of the solution that always includes Nod + 
meniscus survives. It follows from the statement made in section [3~4l that spheres Sph+, n > 0, cannot 
exist when d < 0, so that only the trajectory that contains Sph^ can be observed in {tp, a} plane. 

For d = (the solid sphere on the plane) we show in Appendix IC.2I that the curvature of the Nod~ is 
negative and diverges as H ~ ip~ 2 as tp — > 0, confirming the result reported in [5]. In case 9\ + 9\ > it 
one can observe similar divergence of positive curvature for Nod + at small ijj. The same Nod + meniscus 
can be found in special case 6>i + #i = 7r when the curvature diverges as H ~ ip -1 - In two last cases the 
whole trajectory is represented by Nod + meniscus (see Figure Ob)). 
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(a) 



(b) 



Figure 12: (a) Plot a(ip) for 0i = 30°, 8 2 = 80° and d = -0.1, i.e., -0* = 25.84°, and (b) four menisci 
Nod" (black), Und (cyan), Und^ (red) and Und^ (green) for tp = 45°, if> = 60°, xp = 70° and V = 85°, 
respectively. 



Another divergent behavior of the curvature is observed when d < (the solid sphere intersecting the 
plane). In this case the menisci can exist only for tp > ip* = arccos(l + d). Formula (IC.15I ) shows that the 
curvature diverges as H ~ (tp — V^*) -1 as tp — > tp*. Depending on the parameters values this divergence is 
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observed for menisci, 



Nod 



H < 0, 6»i + 6»i + V>* < 7r, and Nod 



> 0, 01 + 6*1 + V* > 7T- 



(7.1) 



When #i + 6*i + '0* = 7r we show in lC.2l that the Nod - meniscus is forbidden while the Nod + meniscus is 
observed in the whole range ip > ip* and its curvature does not diverge. In Figure [l2]the a(i/j) trajectory 
for contact angles 0\ = 30°, 62 = 80° and d = —0.1 is drawn. We present also the four different menisci 
observed in this case. 

The case of intersecting bodies d < does not lead to change of formula (12.101 ) for meniscus curvature 
H and (12.1 11 ) for surface area S. However formula for the volume V- (d) of meniscus with d < reads 



where integral J s is given in (12- 12b - 

8 Concluding Remarks and Open Problems 

Extending the rigorous approach used [5] to describe the menisci shapes between two touching (d = 0) 
axisymmetric solids, sphere and plane, we develop a theory of pendular rings in its general form for the 
separated (d > 0) or intersecting (d < 0) solids. The main results are listed below. 

1 . The YL equation (12.11 ) with boundary conditions can be viewed as a nonlinear eigenvalue problem. Its 
unduloidal solutions exhibit a discrete spectrum and are enumerated by two indices: the number n G 
Z + U {0} of inflection points on the meniscus meridional profile Ai and the index s = sgn(2a n — 1) 
determined by the shape of a segment of the curve M. touching the solid sphere: the shape is either 
convex, s = 1, or concave, s = — 1. 

Menisci shapes z^(r) and their curvatures play a role of eigenfunctions and eigenvalues of equa- 
tion (12.11) . respectively. The Neeman boundary conditions, two contact angles 6\ and 62, and a single 
governing parameter d > —2 together with one more parameter, the filling angle tp, completely 
determine the meniscus shape z^(r) and curvature evolution H = H^(tp). 

2. For the fixed 6\ and 6 2 the set of the functions H^(tp) behaves in such a way that in the plane {rp, H} 
there exists a bounded domain B where H^(tp) do not exist for any d. 

Under non-linear transformation a n (tp) = H^(ip) sintp/ sin(#i + ip) this domain B in the plane 
{ip, a} takes a simple shape with a smooth boundary /3 ± (V') which we call a balloon. At this bound- 
ary and at the line a = 1/2 there occur all transitions between different types of unduloidal menisci. 




d < 



(7.2) 
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The other two lines a = 0, 1 are also the locations of special menisci: all spherical menisci Sph+ at 
a = 1, all spherical menisci Sph~ at a = and also catenoidal menisci Cat at a = 0. 

A behavior of a n {iji) curves reminds in some cases the 2-dim dynamical system with trajectories 
ramified in a non-simply connected domain A'. This global representation allows to classify possible 
trajectories and introduce a saddle point notion into the PR problem. We observe several types of 
saddle points and give their classification. 

3. If the distance between the solids is non-positive, d < 0, then a single (possibly disconnected) 
sequence of solutions (menisci), that always includes the Nod + type, survives. We describe the 
asymptotic behavior of the mean curvature Hq(iJj) of nodoidal meniscus in vicinity of a singular 
point ip* = arccos(l + d); such singular point does not exist for d > 0. 

Beyond the scope of the present paper we have left several questions which are related to the theory devel- 
oped here. Below we mention two of them. 

1. The theory of pendular rings in the special cases of boundary conditions: 

• #i = 0, this implies a = H and asymptotics dC12| ) of the Nod + meniscus curvature fails. This 
is the case of completely wetted sphere. 

• 02 = 7r/2, the balloon B is reduced to a line a = 1/2 with one singular point tp = tt/2 — 0\. 
This is the case of two solid spheres of equal radii. 

• 0\ = 02, the balloon B is located in region < tp < it — 20\ and the domain A' becomes 
simply connected. This is the case of two solids of the same material. 

• 0\ + 02 = ir, O2 > 0\, the balloon B is located in region < ifi < O2 — 0\ and the domain A' 
becomes simply connected. The physical meaning of this case is unclear. 

2. Stability of pendular rings. In this regard, we know only two papers where the stability was studied 
for the Nod + and Und^ menisci [8] and those menisci which occur as volume decreases from a 
convex bridge [9] between two solid spheres of equal radii. The rest of meniscus types including 
those with inflection points are open for analysis. 
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Appendices 



A Menisci Formulae 

A.l Catenoids Cat 

The catenoid case is the simplest one that requires a solution of the equation (12.21 ) with H = 0. This shape 
corresponds to the transition from the concave nodoid to concave unduloid. The solution for positive x and 
C reads x = C / sinl Using the boundary condition we find 

C = sinii sin^ . (A.l) 

Employing (12.41) we obtain the vertical component of the catenoid meridional profile 

y = [ tantdx = -C f = C fin cot - - In cot — \ . (A.2) 
Jt 2 4 smt V 2 2/ 

Using the boundary condition (12.31 ) we find a relation 

(■. ti h\ sin £2(1 + cosii) 

1 + d — cos ip = sin t\ sin w In cot In cot — = sin t\ sin w In , (A. 3) 

\ 2 2 J smti(l + cost2j 

that implicitly defines the value of the filling angle ip at which catenoid is found. The last condition can be 
rewritten in the form 

tan ^ cot ^ = exp ( - 1 + C ° S ) . (A.4) 
2 2 \ sm ip sin t\ J 

Shape of catenoid is found in parametric form 

x(t) = — , y(t) = C (In cot £ - In cot ^ J . (A.5) 
suit \ 2 2 J 

The solid of rotation volume is given by the formula (12.121) with x = C / sin t. From (12.41) we find dy = 
—Cdtj sint, so that the ring volume reads 

3 f h dt vrC 3 / cos tx cos t 2 , , *i . , t 2 
7TU / q = — — — ^ h In cot — — In cot 



l t2 sin 3 i 2 Vsm 2 ti sin 2 t 2 ' 2 2 

The volume of the Cat meniscus reads 

7rC 3 / cos ti cos i2 ti to \ vr . , , 

F = — - — tt- 5-^ + In cot - In cot -=- - -(2 - 3 cos ip + cos 3 ip) • (A.6) 

2 \ sin t\ sin t2 2 2/3 

The surface area of the Cat meniscus is given by (12.111 ) with x = C / sint producing 

~9 Z"* 1 ^ _ / costi cos to U , to\ 

S = -2ttC 2 — — = tiC 2 [ — + In cot --In cot— . (A.7) 

J t2 sin 3 t Vsin 2 ti sin 2 1 2 2 2 J 
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A.2 Nodoids Nod ± 

The curvature expression reads 

2HV = hih, t 2 ) + sh{t 1 M) ■ (A.8) 
The meniscus nodoid shape is given by 

x(t) = ^ (sint + syWi + c) , y{t) = — [h(t,t 2 ) + sl 2 (t,t 2 )] . (A.9) 

Using the formulas (12.121 ) and (12.221 ) we have for the nodoid volume 

V = 8W [ 4J 3(*i^2) + c/i(ti,t 2 ) -scl 2 (h,t 2 ) +4sJ 2 (ti,t 2 )] - ^(2-3cosV + cos 3 V) • (A.10) 

Using the formulas (12.1 II ) and (12.191 ) we find the nodoid surface area 

S = ^K s (h,t 2 ) . (A.ll) 

A.3 Unduloids Undo 

The curvature expression reads 

2H^ = h(h,t 2 ) + sl 2 (h,t 2 ) . (A.12) 
The meniscus unduloid shape is given by (12.51 ) as 

x(t) = — (sint + sV sin 2 1 + c ) , = — t 2 ) + sl 2 (t, t 2 )\ . (A.13) 

Using the formulae ( 12.12b and ( 12.22b we have for the unduloid volume 

V o = ^777^3 i 4 Mh,t 2 ) + ch(h,t 2 ) - scl 2 (h,t 2 ) + 4sJ 2 (t 1 ,t 2 )]-^(2-3cosV + cos 3 V) . (A.14) 

Using the formulae ( 12.11b and ( 12.19b we have for the unduloid surface area 

S ' = 2jW^ Ks{tut2) ■ (A - 15) 
A.4 Inflectional Unduloids Und^ with Single Inflection Point 

The inflection point u* of unduloid satisfies a condition u 2 + c = for negative c. The inflection point at 



the sphere surface corresponds to 

m = ——. (A. 16) 



sinti 

n 2sim/r 



while when this point is at the plane 



sin t\ + s V sin 2 1 1 - sin 2 i 2 

= ^~ : ; • (A.l/) 

2 sm ?/> 
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Substitution of (IA.16I IA. 17b into (IA. 12b for given s generates equations for the critical values tp J and ijj^ of 
the filling angle at which the inflection point is at the sphere and at the plane, respectively. In case = 
inflectional unduloid reduces to the cylinder reached for t 2 = 7r/2 at t* = n/2 for all Q\ < it/2. It has 
curvature equal to Hf = (2 cos 9\)~ l . 

The integrals in (12.61) and (12.101 ) in case of inflectional unduloid should be broken into two integrals. 
The meridional profile is made of two unduloid profiles matching at the point i.e., u = it*. 

Consider the case of the Und^ meniscus when the profiles touching the plane and the sphere have positive 
and negative curvature, respectively. Using (1A.13I) we write for the convex unduloid part 

x(t) = -^—(smt + V^t + ~c), y(t) = — — [h(t,t 2 ) + h(t,t 2 )} , te{t 2 ,U}. (A.18) 
2ii^ >■ / 2H^ 

The upper concave unduloid part is given by 

x ^ = TTtW fsint- y/ B m 2 t + c + A x ) , y{t) = — — [h(t,t 2 ) - I 2 (t,t 2 ) + Ay] , t e {t x ,U}. 
2H^ V / 2H^ 

The values of A x and A y have to be found from the matching conditions at t = £*. Using u* = sint* = 

yf—c we get t~ = ir — arcsin \/~c and cos t~ = —y/l + c. At the inflection point we find 

\/ = c _ \/^c + A x 
X * ~ 2Hi ~ 2H{ ' 

V* = 7 ^[h(t:,t 2 )+I 2 (t-,t 2 )]=-^-[I 1 (t:,t 2 )-I 2 (t-,t 2 ) + A y \. (A.19) 

The matching conditions produce A x = 0, A y = 2I 2 (t~,t 2 ), leading to the following shape of upper 
concave unduloid 

x(t) = (sint - + 2 , y(t) = -^—[h(t,t 2 )-I 2 (t,t 2 ) + 2I 2 (t~,t 2 )]. (A.20) 

2ii^ v / 2H^ 

Using the second equation in (1A.20I) we have for t = t% 

2H{ V = I 1 (t 1 ,t 2 )-I 2 (t l ,t 2 ) + 2I 2 (t~ ,t 2 ). (A.2 1) 

The case of the Und^ meniscus when the profile touching the plate has negative curvature, and the profile 
touching the sphere has positive curvature is treated similarly and we obtain the general expression for the 
curvature: 

2H(* = h(h,t 2 ) + s[I 2 (t 1 ,t 2 )-2I 2 {t°,t 2 )} . (A.22) 

A.4.1 Shape 

The meniscus shape is given by the following general expressions: 

x (t) = 77777 fsint - sVsm 2 t + c) , y(t) = — [h{t,t 2 ) - sl 2 (t,t 2 )} , t G {t 2 ,t s J, 



2H{ V J 1 a w 2H{ 
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x(t) = ^ (sin t + sVsmH + c^j , y(t) = JL [ h ( t , t 2 ) + s(/ 2 (t, t 2 ) - 2I 2 (tt,t 2 ))] , t G {ti,^}. 
A.4.2 Volume 

As the Und7 meniscus is made of two menisci having shape of concave (upper) and convex (lower) undu- 
loids, a solid of rotation volume equals sum of volumes V/ and V u of lower and upper parts, respectively, 

Vi = J-vi [ 4J 3(**,*2) + ch(t:,t 2 ) -ci 2 (t:,t 2 ) + 4j 2 (t-,t 2 )] , 

V u = ^^[4J 3 (t 1 ,t^)+cI 1 (t 1 ,t-) + cI 2 (t 1 ,t:)-4J 2 (t 1 ,t-)]. 
Adding up the above expressions we have for the meniscus volume 

Vf = - ^ - - { 4 J 3 (t i , t 2 ) + eh (t i , t 2 ) - c [I 2 (t- , t 2 ) - 1 2 (t x , t~ )] + 4 [J 2 (C,i 2 )- J 2 (*i,C)]} 
- ^ (2 - 3 cos V + cos 3 tp) . (A.23) 
The general formula for the volume of inflectional unduloid reads 

v{ = - 7 ^{4j 3 (ti,t 2 ) + c/ 1 (t 1 ,t 2 ) + SC [/ 2 (t:,t 2 )-/ 2 (t 1 ,^)] 

- ±s[J 2 (ti,t 2 )- J 2 (t 1 ,^)]}-|(2-3cos^ + cos 3 V) . (A.24) 

A.4.3 Surface Area 

Apply the same approach to calculation of the surface area of the Und7 meniscus. The area £7" equals the 
sum of the surface areas Si and S u of lower (convex) and upper (concave) parts, respectively, 

Adding up the above expressions we have for the meniscus surface area 



7T 

2(flf)= 



S i~ = 77777^ [*+ (** . *2) + K- (*i - **~ )] - (A-25) 



and we find the formula for the inflectional unduloid surface area 



St = ^sy[K^s(tl,t 2 ) + K s (t 1 ,tl)}. (A.26) 
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A.5 Inflectional Unduloids Undf with Two Inflection Points 

In previous section we consider the simplest basic inflectional unduloid structure characterized by a single 
inflection point. We show that if the inflection point t~ originates at the plane the Und^ meniscus emerges, 
while separation of the inflection point i+ from the sphere generates the Und^ meniscus. 

It is shown in section H]that the value a = 1/2 is a critical point at which a transition between the Und^ 
and the Und^ menisci takes place when t~ = t\ and the single inflection point reaches the solid sphere. 
What does happen when t~ = it — t\ and the inflection point is inside the meridional profile of the Und^ 
meniscus? In this case the second inflection point, namely, with i+ = h, separates from the sphere and we 
observe a meniscus Und^ having two inflection points t+ and t~. The profile of such meniscus is made of 
three unduloid segments - two convex (touching both the sphere and the plane) and a concave one between 
them. 

Consider derivation of the equation for the curvature for this meniscus. Using dA.13b we write for the 
lower convex unduloid part touching the plane 

x(t) = —r (sint+ yWt + c) , y(t) = — — [h(t,t 2 ) + I 2 (t,t 2 )] , t e {t 2 ,t~}. 
The middle concave unduloid part is given by 

x tf) = f sint ~ V^in 2 t + c + A xi ) , V(t) = —— [h(t, t 2 ) - I 2 (t, t 2 ) + A yl ] , t e {t~ , t+}. 

Finally, for the upper convex unduloid part touching the sphere we write 



(*) = 7777+ f sint + Vsm 2 t + c + A x2 ) , y(t) = —— [h(t, t 2 ) + I 2 (t, t 2 ) + A y2 ] , t G {h, t+}. 



2 

The values of A x i and A y i have to be found from the matching conditions at t = ijjf producing 

A xl = A x2 = , Ayi = 2I 2 (t" , t 2 ) , Ay 2 =2I 2 (t-,t 2 )-2I 2 (tt,t 2 ) = 2I 2 (t-,t+), 

leading to the following shape of upper convex unduloid 

X W = 777J+ f sint + Vsm 2 t + c) , y(t) = —— [h(t,t 2 ) + h(t,t 2 ) + 21 2 {t~ ,t+)] ■ (A.27) 
2H2 ^ ' 2J12 

Using the second equation in (IA.27I) we have for t = t\ 

= h(ti,t 2 ) + I 2 (h,t 2 ) + 2I 2 (t-,tt). (A.28) 

The case of the Undj meniscus when the profiles touching both the plate and the sphere have negative 
curvature is treated similarly to obtain 

2H 2 * = h(h,t 2 ) - I 2 (h,t 2 ) + 2/ 2 (*-,t+). (A.29) 
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Thus, we obtain the general expression for the curvature for the Undj menisci 

2fT|tt = Ii(t 1 ,t 2 ) + 8l2(t 1 ,t2) + 2I 2 (t-,ti), (A.30) 
As shown in |D.2| the value J 2 of integral I 2 (£~ , t+) reads I2 = 2- v /^aE(l + 1/c). 

A.5.1 Shape 

The meniscus shape is given by the following general expressions: 

x ( f ) = 77777 (sint + svW 2 t + c), y(t) = — - [Ji(t,t 2 ) + sl 2 (t,t 2 ) ] , t e {t 2 ,t~ s }, (A.31) 



= 77777 f sint - «vWt + c) , y(t) = — [h(t,t 2 ) - s(I 2 (t,t 2 ) - 2/ 2 (C i 2 ))] , * G 



2fT| V / v ' 2iT. 



X W = 77777 (sint + svWt + c ) , y(t) = — - h(t, t 2 ) + sl 2 (t, t 2 ) + 2I 2 ,*€{*!,£} 



2#f V / v ' 2Iff 

A.5.2 Volume 

Inflectional unduloid Und^~ meniscus is made of three menisci having shape of concave (middle segment) 
and convex (upper and lower segments) unduloids; its volume equals the sum of the volumes Vi, V m and 
V u of lower, middle and upper parts, respectively: 



Vl = or J+v* f 4J 3(^, t 2 ) + cfi(C, t 2 ) - cJ 2 (C, t 2 ) + 4J 2 (t- t 2 )] , 

8(Hr ry 

V ™ = 777777C77 [4J 3 , t~) + cli (t+ , t~ ) + d 2 (t+ , t" ) - 4 J 2 (i+ ,£")], 



8(fr+v»' 



'2 

^ = ^7"FjW t 4J 3 (*i.C) + c/i(*i,t+)-c/ 2 (ti,t+)+ 4 J 2 (ti, t+)] • 

8 (^2 J 

Adding up the above expressions we have for the meniscus volume 

y 2 = ,J +N , {4j 3 (ti,f2) + cfi(ti,t2) - cfcccfe) - h(ti,t-) + h{tl,ti)} 
8(Hr ry 

+ 4[J 2 (t-,t 2 )- J 2 (t+ + J 2 (ti,t+)]}-^(2-3cosV + cos 3 V') • (A.32) 
Using the properties of the integrals J 2 and J 2 we find: 



y 2 + = 77r7^f77{ 4 > / 3(ii,t2) + c/i(ti,t 2 )-c[7 2 (ii,t 2 )-2/ 2 ] 



+ 4[J 2 (ti,t 2 )-2J 2 ]}--(2-3cos^ + cos 3 ^) , (A.33) 
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where the value J 2 of integral J^t^ ,t£) is computed in (ID. 14b . The general expression for the menisci 
volume reads 

V 2 = ^7^{ 4J 3( t 1^2) + chihM) -Sc[I 2 (t 1 ,t 2 ) + 2si 2 } 

+ As[J 2 (t 1 ,t 2 ) + 2sJ 2 ]} -^{2-3cosiP + cos 3 ^) . (A.34) 

A.5.3 Surface Area 

Similar approach is applied for calculation of the surface area of the Und^ meniscus. The area S 2 equals 
the sum of the surface areas Si, S m and S u of lower, middle and upper segments, respectively, 

^ = /rT+\9^+^* '* 2 ) ' = — T77+T9 ** ) ' S u = , + ,r,K+{tl,t*). 

2{H 2 y 2(H 2 Y 2{H 2 Y 

Adding up the above expressions we have for the meniscus surface area 

5 2 + = [K + (t:,t 2 ) + K_(tt,t-) + K + (t 1 ,tt)] . (A.35) 

We give a general expression for the menisci surface area 



sr- 71 



2 2(tf|)2 



(A.36) 



A.6 Inflectional Unduloids Undf with Three Inflection Points 

In IA.5I we consider the inflectional unduloid Und^ with two inflection points for the values of inflection 
points parameter < ff < it. There exist also menisci with larger number of inflection points. 

Consider first an inflectional meniscus Und^ with three inflection points. A natural way to generate it 
is to consider the Und^~ unduloid having two inflection points t^ = t\ and t~ = tt — t\ and allow at 
a = 1/2 a third inflection point t~ t = it — 1\ to separate from the sphere. This point appears due to vertical 
translational periodicity of the meridional profile. The profile of such a meniscus is made of four unduloids 
- two convex (one of them touches the sphere) and two concave (one touches the plane). In all formulas 
below we drop additional indices k of tf k . Derivation of the curvature equation for this meniscus is similar 
to the one of Und^ and Und 2 , 

2H^ = h{h,t 2 ) - I 2 (h,t 2 ) + 2/ 2 + 2I 2 (t~,t 2 ). (A. 37) 

A dual inflectional meniscus Und^ is generated from the Und^ unduloid at a = /3 + by separation of a 
third inflection point i+ = t\ from the plane: 

2H+^> = /i(ii,t 2 ) + h(h,t 2 ) + 2/ 2 - 2/ 2 (t+, t 2 ). (A.38) 
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Merging (1A.37IA.38I) we arrive at the general formula for the curvature 

2#f¥ = h(h,t 2 ) + a/ 2 (ti,t 2 ) - 2sl 2 (tl,t 2 ) + 2I 2 . (A.39) 

A.6.1 Shape 

The shape of the Und| meniscus is given by the following general expressions 



x(t) = — (sint - sVsm 2 t + c) , = t 2 ) - sJ 2 (t, t 2 )], iE{i 2 ,i*}, (A.40) 

(*) = ^7 (sint + s\/sin 2 t + c) , j/(t) = ^[h(t,t 2 ) + sl 2 (t,t 2 ) - 2sl 2 (tt,t 2 )}, t G 
z(t) = ^ (sin* - svWt + c) , i/(t) = ^[h{t,t 2 ) - sl 2 {t,t 2 ) + 2/ 2 ], t G {C,**}, (A.41) 



= f Sm * + S ^ sin2 * + C ) > 



2tf| 

= ^[ I i(*^2) + s/ 2 (i,t2) + 2/ 2 -2 S / 2 (^,t 2 )],t E{ti,t;}. (A.42) 

A.6.2 Volume 

Inflectional unduloid Und| meniscus is made of four menisci having shape of concave and convex undu- 
loids; its volume V3 5 equals the sum of the volumes of corresponding parts: 

V 3 = W7^{^Mh,t 2 ) + cI 1 (t 1 ,t 2 ) + sc[I 2 {t 1 ,t 2 )-2sI 2 -2I 2 {t 1 ,tt)] 

- 4s[J 2 (t 1 ,t 2 )-2sJ 2 -2J 2 (t 1 ,^)]}-|(2-3cosV' + cos 3 V). (A.43) 

A.6.3 Surface Area 

The surface area 5| of the Undt meniscus equals the sum of the surface areas of corresponding segments: 



^3 



7T 



2(^1)2 



K- a (£,t 2 ) + K + -K-+K a (t 1 X) 



(A.44) 



A.7 Inflectional Unduloids Undf k with Even Number of Inflection Points 

Generalization of the menisci Undp and Und| to arbitrary even number of inflection points is straightfor- 
ward, and we present here the final formulas for these menisci. Shape of upper unduloid segment touching 
the sphere for the Und| k meniscus reads 

= ^ETT f sin 1 + s V 7 sm 2 t + c) , y(t) = -^—[h(t,t 2 ) + sl 2 (t,t 2 ) + 2kl 2 ] . (A.45) 
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Using the second equation in (1A.45I) we have for t = t\ 



2H s 2k ® = h(t 1 ,t 2 ) + sl 2 {ti,t 2 ) + 2kl 2 . 



(A.46) 



A.7.1 Shape 

The meniscus shape is given for < n < k — 1 by the following general expressions : 



x(t) 



1 



2H. 



2k 



suit + s\/ sin 2 1 + c) , y(t) 



1 



2H. 



2k 



h(t,t 2 ) + sl 2 {t,t 2 ) + 2nl 2 , t G (i 2 ,t- s j(A.47) 



x(t) 

v(*) 



(sint - sVsin 2 t + cj, t e fa ,t+} 



h(t,t 2 ) - s(I 2 {t,t 2 ) - 2I 2 {tT,t 2 ))+2ni 2 



2H, 



2k 



(A.48) 



x(t) 



where 



2H 2k 



s'mt + sy sin 2 1 + c) , y(t) 



2H 2k 



h(t,t 2 ) + sl 2 (t,t 2 ) + 2kl 2 , te{t u tt}{AA9) 



h = t 2 5 n0 + tl{l - 5 n0 ), 



5ij denotes the Kronecker delta. 



A.7.2 Volume 



Volume of inflectional unduloid Und| k meniscus is computed as 



V. 



7T 



2A- 



8(^.) 3 



{4J 3 (ti,t 2 ) + ch(h,t 2 ) - sc[I 2 (h,t 2 ) - 2ksl 2 ] 



7T 



4s[J 2 (h,t 2 ) + 2ksJ 2 }} - ^■(2-3cos^ + cos 3 V') , 



(A.50) 



A.7.3 Surface Area 



Surface area of the Und| k meniscus reads 



5. 



7T 



2k 



K s (t: s ,t 2 ) + k(K + -K-) + K s (t!,t: s ) 



(A.51) 



A. 8 Inflectional Unduloids Und2 k+1 with Odd Number of Inflection Points 

Generalization of the menisci Und^ and Und| to arbitrary odd number of inflection points is straightforward, 
and we present here the final formulas for these menisci. We find for the curvature of the Und| k+1 meniscus 



2H s 2k+1 ^> = h(h,t 2 ) + sl 2 (t h t 2 ) - 2sl 2 {tl,t 2 ) + 2kl 2 . 



(A.52) 
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A.8.1 Shape 

The shape of the Und| k+1 meniscus is given by the following general expressions for < n < k 

1 



where 



x(t) 

y(t) 



2H; 



2k+l 
1 



sin 



t - s\/sin 2 t + c) , te{t 2 ,tl} 



2H. 



[h(t,t 2 )-sl 2 (t,t 2 ) + 2nl 2 



2k+l 



x(t) 



y(t) 



2H, 



I (sint + sVs'm 2 t + c \ , t e {ii,t$,} , 



2k+l 
1 



9 TJ S 
Ln 2k+\ 



h{t,t2) + sl2(t,t2) + 2nl2-2sl2{tt,t 2 ) 



h = hSnO + t* "(I - $no), h = tlSnk + t* S (l ~ S, 



nk ) • 



(A.53) 



(A.54) 



A.8.2 Volume 



Volume of inflectional unduloid Und| k+1 meniscus reads 



V, 



7T 



2fc+l 



,,, Ir , ^{4J3{h,t 2 ) + ch(ti,t 2 ) + ac[J 2 (ti,t 2 ) - 2fcs/ 2 - 2/ 2 (ti,t|)] 



7T 



4s[J 2 (ii,t 2 )-2fcsJ 2 -2J 2 (t 1) ^)]} - ^-(2-3cosV + cos 3 V') • 



A.8.3 Surface Area 

Surface area of inflectional unduloid Und| k+1 meniscus is computed as 

K- S (tt ,t 2 ) + k(k + -k_) + K s (h , t$ ) 



(A.55) 



7T 



°2fc+l — (tts \2 

A.9 Unduloid General Formulas 



(A.56) 



Merging the expressions (1A.46I lA32l) we write the general expression for the curvature of Und^ unduloid 

1 — cos irn 



2H S J> = h(h, t 2 ) + sl 2 (h,t 2 ) + nl 2 - s- 



[h(tt ,t 2 ) + hit' ,t 2 )] ■ (A.57) 



It can be checked by direct computation that the expression in the square brackets in (1A.57I ) evaluates to 
2I 2 (7r/2,t 2 ), and we have 



2H^S> = h(h,t 2 ) + sl 2 (h,t 2 ) + nl 2 - s(l - cos7rn)/ 2 (7r/2,t 2 ), 



(A.58) 
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Replacing in (1A.55I) k by (n — l)/2 we find for the expression in curly brackets 

4J3(ti,*2) + clx(ti,t 2 ) + scl 2 (ti,t 2 ) -ncl 2 + 2scl 2 (ir/2,t 1 ) - 4sJ 2 (tx,t 2 ) + 4nJ 2 - 8sJ 2 (n/2,t 1 ). 
Combining it with (IA.50I ) we find 



Vn = g7^3 i 4 Mh + cli{ti,t 2 ) + £nJ 2 - a{n)[cl 2 (ti,t 2 ) -4 J 2 (ti,t 2 )] + ncl 2 cos ?rn 



+ s(l-cos7rn)[c/2(7r/2,ti)-4J 2 (7r/2,ti)]}--(2-3cosV + cos 3 ^) , (A.59) 



where <r(n) = s • cos 7m. From (IA.5 1 b and (|A.56b one finds 



S ,, 



7T 



K a{n ) U y ' ,h) + [n/2\ (K+ - K-) + K s [t u t 



(A.60) 



where [x\ denotes the floor function. 

B Spheres Sph^ 

In the classical menisci sequence the transition from convex unduloid to convex nodoid takes place through 
formation of a spherical surface Sph^. The inflectional unduloids Und+, n > 0, in the limit a n — > 1 
transform into the surfaces Sph+ made of several spherical segments. The same time the unduloid menisci 
Und~ at small filling angles approach another type of spherical menisci Sph^. Below we treat them both 
as a limiting case of corresponding Und^ menisci. 

B.l Asymptotic Behavior of U nd^ Menisci in Vicinity of Sph^ 

Consider the relation (IA.58I ) and find dependence H = H(tp) for the Und„ meniscus in vicinity of c = 
that can be reached for either a = (for s = — 1, ifi = (f>~ = 0) or a = 1 (for s = 1, ip = Using (12.71 ) 
find asymptotic for c(ip) 

c(V) =c'(^)(V-C) , c / (O = 4 S |2a n -l|a / n (C)sin 2 (0 1 + ^) , sgn c'«) = s, (B.l) 
where > 0. 



To find the asymptotics of the general terms ^fc E(t, \J—l/c) and y/c F(t, \J—l/c) used in (IA.21I) 
we make use of asymptotic expansions for the elliptic integrals based on relations found at iflOl and obtain 

for t < it, 

y/cE(t, y/-l/c) C ~° E(t,c) = 1 - cost - | ^lncot 2 l - - 41n2 + In |c| - 1^ , (B.2) 
y/c~F(t,y/-\/c) C ~°F(t,c) = -| ^lncot 2 ^ -41n2 + ln|c|^ . (B.3) 
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Using the above expressions we find an approximation 

h{tiM) = h{hM)-^M, M = lnftan^cot|Y (B.4) 
The last relation leads to 

/2 = I 2 (tr,i*) = 2-cln2 + |lii(-c), / 2 (7r/2,i 2 )=cost 2 -|mfcot|Y (B.5) 
Substitution of (1B.4IIB.5I ) into (IA.58I ) produces 

on 

2H^((f) s n ) = 2n- (1 + s) cos(0i + <j) s n ) - (1 + s cos vrn) cos 2 + — ln(-c). (B.6) 

The curvature at the sphere is found as 

2n - (1 + s) cos(fli + - (1 + s cos 7rn) cos 2 
ffn(0n) = 2f(^) • (B - 7) 

It follows from (IB.7I) that for the Sph~ sphere the curvature is independent of 9% and reads 

ff-(0) = 2n - (1 - C 2 ° S7rn)c ° sg2 . (B.8) 
For s = 1 we find a condition on the angle 0+ at which sphere Sph+ is observed 

2(1 + d — cos sin(0i + 0+) = sin </>^[2ro — 2 cos(0i + </>^) — (1 + cos Tin) cos 2 ], 
which reduces to 

2(1 + d) sin(0i + 0+) - 2sin0i - sin</>+[2n- (1 + cos vrn) cos 2 ] = 0, (B.9) 
Differentiating the relation (IB .61 ) we obtain in the leading order 

^fr- =i5§T[- 1 + ta <- c )i' (B10) 

Combining ( IB. 10b with the last formula in dB.ll) we find sgn (dH^/dip) = — s at if) = 4> s n . General expres- 
sion for derivative of a n reads 



c«/> sm ti sin ti 



Using it with (IB.8I ) we find 



ff n (0) 2n — (1 — cos 7rn) cos t? 2 



«n = £ n = ■ - = , ■ , , (B.12) 

sm 0i 2a sm 0i 
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implying that for small ip the slope £ n of a increases with growth of index n and decreases with growth of 
the distance d. Applying (1B.2IB.3I ) to integral (12.161) we find its asymptotics ^3 (ti , £2) = cM and obtain 

K s (ti,t 2 ) = 2(1 + s)(cost 2 - costi). 

Thus the expression in square brackets in the general formula for unduloid surface area (1A.60I ) is indepen- 
dent of c and reads 

K = 4n - 2(1 + s cos irn) cos 9 2 - 2(1 + s) cos(6>i + <j) s n ) = AH^((p s n )^((f) s n ). 

Thus the surface area reads in the leading logarithmic order 

*K 27rttQ«) 

and we find that sgn (dS^(ip) / dip) = s at ip = cp s n . The explicit expression for the surface area of Sph„ 
reads 

os / is \ _ Vrn) ,t> ta\ 

nWnJ 2n-(l + a)cos(0 1 + ^)-(l + acos7rn)cose 2 ' ' 

and we find for Sph~ 

Sn(0) = - - a —- -=-. (B.15) 

In — (1 — cos Tin) cos v 2 

Turning to computation of the menisci volume asymptotics we first use (1B.2I IB.3b in integral (12.211) to find 

Mh,t 2 ) = J 3 (ti,t 2 ), J 2 = 4/3. (B.16) 

Using this relation we obtain for the expression V in curly brackets in dA.591 ) 

4 

V = - [4n + 3J 3 (ti, t 2 ) + 3s J 3 (ti, tt/2) + 3s J 3 (tt/2, t 2 ) cos irn] , 
which is independent of c. Thus the volume reads in the leading logarithmic order 

< = W£? ~ v " = ~ ? (2 " 3co8 * + oos3 *> ' (B ' 17) 

and we immediately find that sgn (dV^ (ip) / dip) = s atip = 4> s n . The explicit expression for the volume of 
Sph„ reads 

4vr^ 3 (<&) [4n + 3 J 3 (t 1 , t 2 ) + 3s J 3 (t 1 , tt/2) + 3s J 3 (tt/2, t 2 ) cos 7m] 
3|2n — (1 + s) cos ii — (1 + s cos nn) cos p 2 J 

with ti = 0i + <p s n and t 2 = it — 9 2 . Recalling that 

J 3 (ti,t 2 ) = [G(t 2 )-G(ti)}/3 , where G(t) = 3 cost - cos 3 t , G(tt/2) = , 
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we have 

v s U s , _ 4vr^ 3 (^)[4n - (1 + s)G(0 1 + - (1 + SCOS 7m)G(fl 2 )] _ 

3[2n-(l + a)cos(6>i + &)-(l + «cos7rn)a*0 2 ] 3 1 nJ ' 

and we find that for Sph~ 

V-frt 4^ 3 [4n-(l-cos^n)G(g 2 )] 
3[2n — (1 — cos im) cos f)2\ 
the volume does not depend on 9\. Using (IB.ll ) and (IB. 101 ) in (IB- lib we find 

, . , , , ^ sin 0i 

= m -F^TU ^ • ■,/,+ - , • (B " 21) 

w — nsm0n smti ln(— cj sm0„smri 

which produces two important formulas for n = 

<W) = ■ » (B.22) 

sin^o sm(6'i + </>q ) 

and for n>l, 

, / <+s sin gi g n 1 

«nWJ - gin 2 ^+ sin 2 ( ^ + ^ ' n l n (_ c ) ~ dsin 3^ l n (_ c ) " 

The last equality in (1B.23I ) makes use of (13.211) . i.e., n</>+ ~ d sin 6\, and ^ ~ d when n — > 00. 
B.2 Sphere Sph^ 

The curvature of Sph~, n > menisci is given by ( IB. 81) . Using the general formula (lA.47l - lA.49b for the 
unduloid Und^ k we find in the limit ip — > that the meniscus is presented by a sequence of coaxial spheres 
of the radius T2k = V-^2fc = ^/ (^) w * tn tne centers located at {0, (2? + l)r2fc} for i = 0, 1, . . . , k — 1. 

Similarly, for odd number n = 2k + 1 of inflection points we have k full spheres of the radius r 2 k+\ = 
l/H^ k+l = dj (2k + 1 — cos 62) and a spherical cap on the plane with the same radius. The centers of the 
full spheres are at {0, d — (2i + l)r2fc+i} for % = 0, 1, . . . , A; — 1 and the center of the spherical cap is at 
{0, -r 2 k+i cos 6 2 }. 

B.3 Sphere Sph+ 

The curvature of Sph+ menisci is given by (IB.7I) with ip = In case n = 2k the meniscus is represented 
by the segments of the spherical surface touching each other at the vertical axis at the points with zero 
abscissa and ordinates (2i + 1 + cos t2)/H£, i = 0, 1, . . . , k. 

In case n = 2k + 1 the meniscus is made of the segments of the spherical surface touching each other at 
the vertical axis at the points with zero abscissa and ordinates 2i/H+. It follows from dB.91 ) that the value 
of the filling angle 4>2k+i ^ oes not depend on the value of the angle 62- 
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In a particular case of wetting sphere 9\ = one finds from dB.9b that Sph+ meniscus exists at tp = tt 
that gives for the curvature = [n + 1 + (1 — cos 7m)) cos t^\/{2 + d). Taking into account that 
= a n < 1 we find an existence condition for this meniscus 

d > n — 1 + (1 — cos 7m) cos £2- (B.24) 

The condition (IB.9b at #1 = is also satisfied for all < ip < tt when d = n — 1 + (1 — cos 7m) cos £2- 

C Special Properties of Nodoids Nod ± 

In this appendix we discuss the extremal properties of nodoidal menisci (local minimum of curvature and 
local maxima of the surface area and volume) and the asymptotic behavior of the nodoidal curvature in 
vicinity of singular point ip* = arccos(l + d). 

C.l Non-monotonic Behavior of Nod + Meniscus Characteristics 

Show that the Nod + meniscus always has a local minimum of curvature and local maxima of the surface 
area and volume. First, we show that the curvature always grows when the filling angle reaches tt, so that 
the derivative dH/difi at ip = tt is positive. In this range the convex nodoid is observed, so that we start 
with the asymptotics of the general terms yfc E(t, y/—l/c) and y/c F(t, \J — 1/c) in (IB.2IIB.3b valid for 
t < 7r and also find for t > tt 

\/cE(t,y/-l/c) C ~° -E{t,c) + 4-c(-l -41n2 + ln|c|) , (C.l) 
y/cF(t,y/-l/c) C ~°-F(i,c) -2c(-41n2 + ln|c|). (C.2) 

We use (IA.8I ) with s = 1 corresponding to convex nodoid and substitute into it the expressions (1B.2I 
IB.3I ) for t = t2 = tt — 62 < tt and (IC.1IIC.2I ) for t = t\ = tt + 6\ > tt. Retaining the leading terms only 
we arrive at 

2(d + 2)#«2(l-cos6> 2 ) + - f 2-81n2 + 21nc + mtan 2 ytan 2 y J , (C.3) 

where c w AH(tt — tp) sin 9\ is positive for ip < tt. As c'(V') < the leading term in the above expression 
is cine which for positive c guarantees curvature growth in the vicinity of ip = tt. As H(0) > H(tt) the 
curvature dependence on the filling angle cannot be monotonous one, and the curvature should have a local 
maximum and a local minimum. Thus, as the curvature at the spherical meniscus Sph^ is a decreasing 
function of the filling angle (see Appendix lB.il ) and the same time in the vicinity of ip = tt it always grows, 
it always has a local minimum at the convex nodoid meniscus. 
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Volume behavior analysis gives in the leading order 

v(! \ - 647r(ri + 2) 3 (2-cos£ 2 )(l + cost 2 ) 2 4tt 7r(d + 2) 3 (2 - cos t 2 ) 4tt 

W 3[4(1 + cost 2 ) + cine] 3 3' W 3(l + cost 2 ) 3' ( ' 

and its derivative in ip reads in the leading order 

dV(ip) __ 7r(d + 2) 3 (2-cost 2 )lnc , 

~ 4(l + cost 2 ) 2 C W ' 

and the volume decreases at ip = n. Comparing the volume at two extreme values of the filling angle we 
find that V(0) < V(ir), which implies that its behavior is non-monotonic and it should have at least one 
local maximum and one local minimum. As the volume grows at the spherical meniscus (see lB.ll ) its local 
maximum is observed on convex nodoid meniscus. 

Analysis of the surface area of the convex nodoid in the vicinity of ip = 7r is done similarly to that of 
performed at small filling angles in Appendix IB. II The area is given by (IA. lib with s = 1 that gives in the 
leading logarithmic order 

= 3Mrf+ 2)»(1 + cos «,) sw = Md + lf (C5) 

yY> [4(1 + cost 2 ) + clnc] 2 KJ l + cost 2 

and its derivative in ip reads in the leading order 

dSjjj) 7r( ( i + 2) 2 lnc 
# ~ (l + cost 2 ) 2 CW " 

Noting that c'(t/j) < we find that the derivative of the surface area w.r.t. the filling angle is negative for 
ip = 7r and the surface area decreases. Comparing the surface area at two extreme values of the filling angle 
we find that 5(0) < S(tt), which implies that its behavior is non-monotonic and it should have at least one 
local maximum and one local minimum. As the surface area grows at the spherical meniscus (see IB - 1 b its 
local maximum is observed on convex nodoid meniscus. 

C.2 Asymptotics of Nodoid Curvature 

Here we discuss divergence of the curvatures of convex and concave nodoids and corresponding asymp- 
totics. Consider equation (12.101 ) for the Nod s menisci, 

2HV = h + sl 2 , (C.6) 

where I\ and J 2 are defined in (12.101 ). (12.131 ) and (12- 14b - For both nodoids c > 0, therefore integrals 
l\ and I 2 are always convergent and divergence appears only when ^ is vanishing. This happens when 
— 2 < d < 0, i.e., the divergence does not occur when the solid bodies are separated. 
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Consider first the case d = for which a singular point is -0^ = and in its vicinity ip <C 1 we obtain 
\£ ~ ^ 2 /2- Choose the power law of divergence, H ~ Uip~^°, where U, A) > 0, then in accordance with 
(12771) we find 



-4ffsini/!siri£i ~ U^ 1 '? , A) < 1, f7i = -4E7'sin0i > 0, 
- < 4Fsin^(FsinV-sinti) = C/ 2 , A) = 1» ^2 = 4f7(Z7 - sin 6>i) > 0, 



(C.7) 



AH 



2 sin 2 



! V =± c/ 3 V' 2(1 ~ /3o) . A) > 1. ^3 = 4f/ 2 > 0, 



that for ip — > implies c ~ 0, A) < 1, c = 4U(U — sin#i), A> = 1 an d c ~ oo, A) > 1. 

Consider the integrals I\ and /2- For the first of them we have I\ ~ J* + ^sinfli, where 7* = 
— (cos #i + cos #2)- Regarding I2, denote P(t, c) = sin 2 t/\/c + sin 2 i and obtain 

_ j;% 2 p(t, u 3 ^-M)dt + i>p(6 x , u^-M) , ^ > i- 

Substituting (IC.8I ) into dC6l ) and making use of asymptotics dB.41 ) we find 

IiJ - ^V 1 "* In (tan f tan f ) + , A) < 1, 

fa P(t, U 2 )dt + ^P(9 1 , U 2 ) , A) = 1, < C9 ) 
vfc^ 1 L% 2 sin 2 i (it + ^P(0i, U^-*>)) , A) > 1, 

Solve equation (1C.9I ) in two cases. First, if #1 + #2 / tt, then preserving the leading terms in tp we get 
for the Nod+ (6\ + 9 2 > vr) solutions, 

u^ 2 -^ = 211 A) < i; ^ 2 ~ /3 ° = A* + / P(t, u 2 )dt, A) = i; ^ 2_A) = ij, A) > 1, 

which yields J7 = I*, (3q = 2. In the case of the Nod~ nodoid and 9\ + 62 < it we get 

J^i-A) In (tan \ tan , A) < 1 

it-lHU.PitMdt, A) = i » ( C1 °) 

A* , /3 > 1 

satisfied for C7 = I*, A) = 2 only. Thus, in the generic setup #1 + 02 7^ n the both nodoidal menisci have 
divergent curvature, 



Uib 



2-00 



H 



COS 9\ + COS 02 

^2 



(C.ll) 



In case 62 = it/ 2 its expression coincides with estimate (11.21 ) derived by simple considerations. 
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Consider a special case 6\ + 62 = vr and 6\ 7^ 0, 7r, for Nod + rewriting dC9l ) in leading terms, 

= 2V>sm0a, /3 < 1; Utf~to= ™* m6 \ ^, ft) = 1; [Z^ 2 "* = ^ sin X , ft, > 1 , 

2t7 — sm#i 

which is satisfied for U = § sin 0i, ft) = 1, i.e., 

ff = ^, flx^O.TT. (C.12) 

In case Nod - with #1 + 62 = vr, 6\ / 0, 7r, we have from (1C.9I ) after substitution of £/i, t/2 and U3 defined 
indCTT) 

^ 2_/3o = -2C/V' 2 " /30 , ft) < 1; Z7V 2 ~ ^ = 2 ~~— Sm °l V sin 6 1 , ft, = 1; Uifj 2 "? =ij, sin0i, ft, > 1 

2U — sm 0i 

The first and third equations cannot be satisfied due to restrictions on ft, and U / 0. The second equation 

U(2U-sm6 1 ) = 2(£/-sin0i)sin0i , ft) = 1 , (C.13) 

does not admit real solutions, so the Nod - meniscus is forbidden in the special case 61 + 62 = vr. 

For — 2 < d < a singular point ip* = arccos(l + d) > does exist and in its vicinity ip — ip* = 
r\ <C 1 we obtain ^ ~ 77 sin i/j* . Choosing the power law of divergence, we find 

c ~ 4y 2 rj~ 2 ^ 1 sin 2 ?/>*. Write the leading in r\ terms of integrals I\ and /2 

/"** 

ii = J** +r7sint*, J 2 = / P(t, AV 2 sin 2 ^rf m )dt + 7?P(^, 4F 2 sin 2 

Jt 2 

where /** = -(cos t| + cos 2 ) and t* = 9i + ip*, and substitute them into (IC.6b . 



2T/r ? 1 - /31 sin V* ^ IT* + ^ sin t\ + S7? - [ / * sin 2 tdt + r) sin 2 t\ | . (C. 14) 

Wsmij)* \J t2 J 

In general case, 6\ + #2 + ip* 7^ ^ we nave f° r b° tri nodoids Nod + (#1 + 02 + V 1 * > vr) and Nod - (0i + 

6*2 + *0* < vr), 

cos (6*i + VQ + cos0 2 1 „ - 

2 'V-V*' 

In case 02 = vr/2 its expression coincides with estimate (11.41 ) derived by simple considerations. 
The special case 6\ + 02 + ip* = vr leads to 

s sin 2 i* 

2y r? -Pi sin^* ~ sintt + — -r/\ 

21/ sm 

satisfied by /3i = and H does not diverge in vicinity of the critical value tp* . This conclusion holds for 
any other (non power law) divergence H ~ g(rj) when 6\ + 02 + ip* = vr, 

s sin 2 t\ 1 



2Vg(rj) sin ^* — sin + 



2V sin-(/>* 
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Show that for 6\ + 62 + ip* = ft the Nod meniscus is forbidden while the Nod + meniscus is allowed 
for ip > ip*. First use (IB.9I) for n = to find the value <p^ at which the sphere Sph^ is observed: 

(1 + d) sin(0i + 4>q) + smtpQ cos 02 = sin0i. 

Direct computation shows that <Pq = ip* = ft — 6\ — O2 satisfies the above equation, so that the sphere 
Sphjj" exists at ip = ip*. As the meniscus Nod - exists in the range ip < ip* = ip* which is forbidden due to 
intersection, we conclude that Nod - cannot be observed in this special case. The same time, the meniscus 
Nod + is allowed for ip > 4>q = ip*. The value of the Nod + curvature at ip = ip* can be obtained by noting 
that it is equal to the Sph^ curvature that reads H = sin(0i + ip*) / sin ip*. The Table 3 (where 77 = ip — ip* 
and ip* = arccos(l + d)) summarizes the asymptotic behavior of the Nod 1 * 1 menisci curvature. 



Table 3. 





d 


6*1 + 2 + ip* < 7T 


61 + 2 + V>* = ft 


01 + 2 + Ip* > ft 


Nod - 


= 


-(cos (0i + ip*) + cos 2 )ri~ 2 


forbidden 


forbidden 


Nod- 


< 


— (cos (0i + ip*) + cos 9 2 )ri^ 1 


forbidden 


forbidden 


Nod+ 


= 




(3/2sin0i)r ? - 1 


-(cos (0i + ip*) + cos 2 )r?" 2 


Nod+ 


< 




sin(0i + ip*)/ sinip* 


— (cos (0i + ip*) + cos 6 2 )ri~ 1 



The empty entries in Table 3 indicate that the Nod + meniscus does not exist in the vicinity of ip* contrary 
to the "forbidden" entry that means that the corresponding meniscus does not exist in the whole range 

tp* < ip < ft- 

D Computation of Elliptic Integrals 

In this appendix we derive formulas for computation of the elliptic integrals used in the main text. 

D.l Conjugation of Elliptic Integrals 

Here we prove that 

f z dt 1 ~T Z Jt 

/ / 5 = ~r / / . (D.l) 

J Vsin z t + c v c Jo 1 _|_ sal* 
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where A(z) stands for complex conjugation of the function A(z). The case c > is trivial and the operation 
A(z) can be omitted there. Consider negative c and rewrite the l.h.s. of (ID. II) as follows 

r z dt 

R(u;0,z) = / c=-u 2 , (D.2) 



o \j sin 2 1 - v 2 



and focus on two cases: 



1. sin 2 t - v 2 < 0, when < t < z, 



2. sin 2 t — v 2 < 0, when < t < z±, and sin 2 1 — v 2 > 0, when z* < t < z. 

In the first case the integral in (ID .21 ) is purely imaginary, 

f z dt 1 /"* 

R(u;0,z)= = — / - (D.3) 

•io Vsin 2 t - v 2 w Jo V 1 - v~ 2 sin 2 t 

where an integral in the r.h.s. of (ID.3I) is real (positive). Thus, equality (ID. 1 1) holds also in this case. In the 
second case write R(u; 0, z) as a sum R(v; 0, z) = 0, z*) + 2*, z), 

f 2 * (it f 2 (it 

R(v;0,z*)= - , R(u;z„z)= - (DA) 

Jo y sin t — v 1 Jz* Ysin t — v 2 

The first integral in dD.41 ) for sin 2 £ — i/ 2 < 0, is purely imaginary and can be calculated using dD.31 ) 

1 /" z * dt 

R(u;Q, Zir ) = — -., sfc = Mv, (D.5) 

±^ ./o Yi - sin 2 t 

Equality (ID. II ) holds for 0, z*). The second integral R(y\ z*, z), where sin 2 t — i/ 2 > 0, is positive, so 
-R(V; 0, z) can be represented as follows, 

dt 1 ( r* dt , . r z dt \ m „ 

± / - . (D.6) 







\/ sin 2 t — v 2 ^ v \Jo \/l — v~ 2 sin 2 t J z+ y 7 sin 2 t — v 2 
Consider now another integral, 

1 f z dt 1 f z dt 

T(v;z*,z) = — - R {y.Q, z )= / (D.7) 

± w Jo Yi - ^~ 2 sin 2 t ±»" Yi - sin 2 t 

which can be rewritten as follows 

1 f z dt f z dt 

T{y-z^z) = —\ -- = - \ ; n (D.8) 



\J v 2 sin 2 1 — 1 A* y / sin 2 1 
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which is a negative number. Comparing the latter with (ID.4I) we obtain T(y\ z*, z) = —R{v; z*, z). 

Combining the last equality with (ID.7I ) we obtain 

1 f z dt if r* fit /■* (it \ 

Tiv ; = . (D.9) 



±^ •/ Vl - z/- 2 sin 2 t \J0 yfl - u~ 2 sin 2 t A* y 7 sin 2 t - v 2 

By comparison (ID .61 ) and dD.91 ) we find finally 



r 

.7 



* 1 r * (D.tO, 



'0 \/ sin 2 t — i^ 2 JO a/1 — v~ 2 sin 2 t 

Keeping in mind that we have taken y/c = ±iv in (ID.5I — lD77l > and (ID.9I — IdTTQT ). we arrive at (ID. lb . 
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D.2 Computation of Elliptic Integrals at Special Limit Values tf 

The integrals I2 = h(t*,t+), I3 = h(t~ , t+) and J 2 = J2(t~,t£) enter numerous formulas for unduloids 
Und„ so that it is instructive to find their explicit expression through the complete elliptic integrals of the 
first K and second E kind. In derivation we used relations from lfTTl[T2l . We start with the general relations 

E{tt + im, k) = k[E{-c) - (c + l)K(-c)} + 2nE{k 2 ), F(t+ + im, k) = ^/^K(-c) + 2nK(k 2 ), 

(D.ll) 

using them with n = 0, — 1 for i+, t~, respectively. From definition (12.151) of I2 integral we obtain for I2 

I 2 = 2^~c[E{k 2 ) - K{k 2 ) - k(E(-c) - K(-c))}, 
where the expression in the square brackets simplifies to iE(l + 1/c) leading to 

J 2 = 2 v /= c£(l + 1/c). (D.12) 

We also find 



(D.13) 



dh cE(l + l/c)+K(l + l/c) 

dc ~ (l + c )y^ 

Using (12.211 ) it is easy to check by direct computation that 

t l + c f E(t~,k)- E(t+,k) 
J2-—I2 + § , 

and we find 

■h = 2(1 » +C) \/ Z ^^(1 + 1/c) + \[%E(-c) - i(l + c)K(-c) - V^E(-l/c)}. (D.14) 
Finally, using (12.161) we find 

h = -2[icK(-c) + ^K(-l/c)}, 

and using lfT2l we arrive at 

h = -2y/=cK(l + 1/c). (D.15) 
Collecting the expressions (ID.12ID.15I) and using the definition ( 12. 19b we find 



K s = K T (t-,tt) = 4VTT^ + 2sV^[2E(l + 1/c) - K(l + 1/c)]. (D.16) 

Finally, consider integral I4 = I^(t~,tf) = I± c + J^, which is written as a sum of a constant term 
I^c and a divergent part Z^. This representation follows from (12.241) where the second term diverges as 
c + sin 2 t^ = and we find 

~ sin 2t+ 2tant+ 



L 4d 



(l + c)Vc + sin 2 tf vc + sin 2 t* 
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Introducing c = (4e 2 — 1) sin 2 t+ where e — > we obtain 

hd = ^-+=^=- (D.17) 

ecosr* ey 1 + c 

Turning to the constant term in I4 we compute the first term in ( 12.241 ) using the relations ( ID. Ill ) and find 

h c = {^(A; 2 ) - j^- c [E(k 2 ) ~ kEi-c)}} • (D.18) 

Using the relation from lfT2l 

E{z) = \[zE(\jz) - iE(l -z) + izK{\ - z) + (1 - z)K(z), 
we find for z = — 1/c 

1 + c 



S(-l/c) - y/-l/cE(-c) = -iE{\ + 1/c) - (i/c)A-(l + 1/c) + K(-l/c). 

c 

Substituting it in (ID. 18b and comparing the result with (ID. 13b we find 

h = -2l' 2 (c) + (D.19) 
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